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Within  the  configuration  interaction  (CI)  framework,  a  higher  order 
perturbation  approach  has  been  developed  v;hich  allows  a  drastic  reduction 
in  the  number  of  n-particle  functions  required  to  give  the  energy  of  a 
CI  wavefunction  to  a  desired  numerical  accuracy.   By  means  of  this  method, 
significant  improvements  over  the  standard  Hartree-Fock  (HF)  and  truncated 
CI  calculations  are  obtainable.   Since  the  procedure  does  not  require 
solutions  of  large  secular  equations,  the  technique  also  has  implications 
for  full  CI  calculations.   As  this  procedure  complements  the  natural  orbi- 
tal methods,  the  two  techniques  can  be  used  in  conjunction  for  additional 
flexibility. 

The  theory  derived  to  effect  the  reduction  process,  i.e.,  the  reduced 
partitioning  procedure,  employs  the  techniques  of  inner  and  outer  projec- 
tions and  Pade'  approximants  to  yield  a  tractable  computational  scheme  with 
kno'.-m  convergence  properties.   The  theory  is  formulated  both  for  a  single 
unperturbed  reference  function  and  for  a  multi-dimensional  reference  space 
to  use  in  applications  to  excited  states.   In  the  latter  case,  a  method 

viii 


which  uniformly  treats  the  excited  states  as  well  as  one  which  permits  a 
selective  reduction  for  one  specific  state  is  obtained.  These  reduction 
equations  are  also  shovm  to  be  related  to  the  method  of  steepest  descent. 

The  first  order  solution  is  of  particular  significance.   In  the 
examples  studied,  this  result,  which  corresponds  to  the  application  of  a 
geometric  sumrule  to  the  perturbation  expansion,  is  found  to  give  sixty 
to  seventy  per  cent  of  the  possible  energy  improvement  obtainable  from 
the  specified  basis.   This  solution  also  has  some  computationally  simpli- 
fying features  when  the  reference  function  is  of  the  HF  or  CI  variety. 

Preliminary  applications  to  the  H  molecule  and  the  HeH^  molecular 
ion  tend  to  justify  the  reduction  procedure  for  more  general  cases.   In 
these  examples,  no  more  than  a  tv/elfth  order  solution  is  required  to  give 
essential  agreement  with  the  full  CI  result  in  each  case,  even  though  as 
many  as  eighty-two  configurations  constitute  the  full  n-particle  basis. 
In  addition,  it  is  found  that  only  a  third  order  solution  is  required  to 
account  for  more  than  ninety  per  cent  of  the  possible  correlation  improve- 
ment over  a  HF  calculation.   These  calculations  presented  here  also  in- 
clude potential  curves  for  the  HF,  the  full  CI,  and  the  perturbation 
corrections  for  selected  states  of  HeH^ , 
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CHAPTER  I.   INTRODUCTION 

The  majority  of  theoretical  studies  of  atoms,  molecules,  and  solids 
has  its  origin  in  the  time  dependent,  non-relativistic  Schrodinger  equa- 
tion [1], 

More  comm.only,  though,  one  focuses  on  its  progeny,  the  time- independent 
form 

HY  =  E  Y   ,  (1.2) 

n    n  n 

which  defines  the  "statioiiary  state"  eigenfunctions ,  (Y  }.   The  complete 
Hamiltonian  operator  for  a  given  problem  is  an  enormously  complex  entity, 
whose  eigenfunctions  would  contain  all  electronic,  vibrational,  and  rota- 
tional motions  that  are  possible  for  the  system  as  well  as  a  variety  of 
hyper-fine  interactions.   This  enigma  is  never  considered  in  its  entirety, 
but,  subject  to  the  problem  of  interest,  one  makes  certain  assumptions  as 
to  the  form  of  H.   Even  then  the  equations  obtained  arc  probably  still 
partially  coupled  and  require  another  approximation  to  effect  a  separation. 
It  is  this  process  and  the  Born-Oppenheimer,  or  "clamped  nuclei  ,"  approxi- 
mation [2]  that  lead  us  to  the  usual  Hamiltonian  Cor  the  electronic  motion, 

H  =  -J2EV^'+  E  1/r..  -  E>:z  /r  .  -1-  S  Z  Z^/R  ^  .     (1.3) 
1     i<j    -■  a  1  Q'<p       "^ 

Tlie  quantities  in  (1.3)  are  the  kinetic  energy  operator  and  terms  arising 
from  the  coulombic  potentials  for  elccti'on-electron,  electron-nuclear,  and 
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nuclear-nuclear  interactions.   Within  the  clamped  nuclcii  appi'oximation 
the  last  term  is  merely  a  constant.   I'/hen  additional  interactions  are  to 
be  considered,  they  are  cvistomarily  treated  as  a  perturbation  of  (1.3). 

Even  though  (1.3)  is  a  distinct  simplification,  it  is  still  impossible 
to  solve  (1.2)  for  a  system  with  more  than  one  electron,  necessitating 
approximate  solutions  as  the  only  recourse.   In  this  respect,  powerful 
methods  have  been  developed  to  approximately  solve  (1.2).   These  fall  pri- 
marily into  two  categories,  variational  methods,  based  on  the  Ritz  varia- 
tional principle  [3],  and  perturbation  theory.   The  former  has  tradition- 
ally been  the  more  important  of  the  two,  with  the  latter  essentially 
limited  to  the  effects  of  external  fields.   However,  a  renev;ed  interest  in 
perturbation  methods  is  very  much  in  evidence  currently. 

Since  (1.3)  is  bounded  from  below,  the  variational  theorem  allows  one 
to  approximate  Y  by  some  Y  which  contains  several  parameters  to  be  opti- 
mized, form  the  functional 

E  =  <y|h|?)/<y|?>  ,  (1.4) 

and  by  minimizing  E  V7ith  respect  to  the  parameters  embedded  in  Y,  determine 
an  approximate  wavefunction  and  an  E  s  E.   With  a  sufficiently  flexible 
expression  for  Y,  one  can  thusly  obtain  as  good  a  solution  as  desired. 
In  this  manner  the  energies  of  a  few  simple  two- and  three-electron  systems 
have  been  calculated  to  experimental  or  better  accuracy. 

For  less  simple  systems,  two  relatively  general  variational  methods 
have  replaced  the  "special  case"  solutions  of  the  previous  period.   These 
are  the  Hartree-Fock  method  [4],  which  for  molecules  is  invariably  imple- 
mented as  formulated  by  Roothaan  [5],  and  the  method  of  superposition  of 
configurations  or  configuration  Interaction.   The  ordinary  Hartree-Fock 
(HF)  wavefvnction  for  an  atom  or  molecule  with  n  electrons  is  assumed  to 


be  a  single  Slater  determinant, 

$^  =/(^-^(l)i'2(2)  ...  ^y^/n))  (1.5) 

composed  of  n  one-particle  functions  (spin-orbitals)  V7here  each  spin- 
orbital  is  a  spatial  orbital  5.  associated  with  either  q-  or  3  spin.   In 
the  Roothaan  procedure,  the  functions  {5.}  ^^rc  assumed  to  be  approximated 
by  expanding  in  a  given  set  of  functions  [X  ]> 

§.  =  E  X  c  .  (1.6) 

[i  ^ 

with  arbitrary  expansion  coefficients  c   .   Tlie  coefficients  can  be 
determined  by  the  variational  principle  (1.4)  to  give  the  lowest  possible 
energy  for  the  single  determinant  form.   Besides  being  conceptually  satis- 
fying, a  wavefunction  of  this  form  possesses  a  degree  of  justification 
for  the  prediction  of  ionization  potentials  and  one-electron  properties 
due  to  the  Koopmans  [7]  and  Moeller-Plesset  [8]  theorems.   On  the  other 
hand,  all  single  determinant  wavefunctions  are  unable  to  properly  treat 
the  instantaneous  electron  interactions,  i.e.,  "correlation,"  which  arises 
from  the  second  term  of  (1.3).   Tliis  has  several  consequences,  one  of  V7hich 
is  that  a  standard  Hartree-Fock  wavefunction  for  a  molecule  cannot  general- 
ly be  used  to  obtain  reasonable  potential  energy  curves  since  it  will  not 
dissociate  properly  except  in  a  few  very  special  cases.   Furthermore,  one 
encounters  difficulties  with  sjiranetry  in  proceeding  to  open  shell  studies 
resulting  in  the  near  impossibility  of  treating  excited  states. 

The  configuration  interaction  (CI)  procedure  effectively  eliminates 
the  complications  inherent  in  the  Hartree-Fock  scheme,  but  not  without 


^  "A  common  basis  set  is  the  Slater  type  orbltals  [6],  X(l,m,n)  = 
r   e   Y(l,m),  v/here  l,m,  and  n  are  the  usual  quantimi  numbers.,  r  the  mag- 
nitude of  the  radial  coordinate  relative  to  an  arbitrary  origin,  Y(l,m)  the 
spherical  harmonics,  and  q'  the  orbital  exponent  which  is  to  be  selected  in 
some  manner. 
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causing  some  new  problems  of  its  ox^m.   If  one  considers  the  expression 

(1.6)  and  if  there  are  m  expansion  functions  [X  },  then  with  spin,  it 

follows  that  2m(iy.)  are  possible.   Consequently,  if  2m  >  n,  many  more 

than  one  determinant  may  be  generated  from  these  {{' .]   by  forming  the 

(   )  combinations.   In  the  CI  method  one  can  consider  these  determinants, 
n 

{d  } ,  to  be  the  basis  for  a  trial  wavefunction, 

■^  =  E  \C^    ,  (1.7) 

k 

with  the  [c  ]  to  be  determined  from  (1.4).   Other  than  the  multi-determinantal 
nature  of  the  CI  approach,  which  permits  electron  correlation,  the  primary 
distinction  between  it  and  Hartree-Fock  theory  is  that  the  expansion  coeffi- 
cients in  the  [§ . }  are  not  usually  determined  via  an  energy  criterion  but 
rather  fixed  --  such  as  by  orthogonalization  of  the  (x  }  --  thereby  placing 
the  flexibility  of  the  wavefunction  exclusively  within  the  [C  } .   (Exclud- 
ing  the  multi-conf igurational  Hartree-Fock  techniques,  the  one  exception  to 
this  statement  involves  natural  orbital  iterations  as  discussed  below.) 
Immediately,  a  problem  with  the  expression  (1.7)  becomes  apparent, 

since  the  number  of  determinants  (   )  rapidly  becomes  astronomical  even 

n 

for  a  relatively  modest  basis  set,  f y  ] .   If  one  is  interested  in  a  state 
or  states  of  a  certain  symmetry,  many  of  these  determinants  can  be  excluded 
by  group  theory  since  they  will  contain  no  component  of  that  state.   Even 
better,  the  individual  determinants  which  do  contribute  to  the  state  of 
interest  can  be  grouped  into  linear  combinations,  {b  },  such  that  each  B 

Li  Xj 

is  a  pure  symmetry  function.   By  so  doing,  the  expansioia  (1.7)  can  be 

substantially  shortened,  but,  even  so,  the  number  of  configurations  is 

still  likely  to  be  very  large. 

If  one  forms  the  approximate  wavefunction 

Y  =  E  B  C   ,  (1.8) 

L  ^  ^ 
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by  including  every  co:if iguration  that  can  be  constructed  from  the  given 
basis  [x  ],  then  the  result  obtained  from  the  variationtil  pi"inciple  is 
termed  the  "full"  CI  solution  for  the  problem.   This  solution  has  some 
convenient  properties.   First,  it  is  invariant  to  any  transformation 
among  the  basis  functions  {x  ]  and  thus  is  independent  of  the  explicit 
form  for  the  {5.};  or,  equivalently,  the  final  result  is  solely  determined 
by  the  initial  basis  [x  }.   The  second  advantage  is  that  by  using  all  pos- 
sible  {b  },  no  dubious  selection  of  configurations  needs  to  be  made.   Un- 
fortunately,  though,  the  number  of  configurations  needed  for  the   full   CI 
is  still  usually  much  too  large  to  enable  one  to  use  all  of  them,  a  fact 
which  leads  to  severe  diagonalization  problems  among  others.   In  order  to 
avoid  most  of  these  problems,  it  is  necessary  to  make  a  somewhat  arbitrary 
selection  of  the  configurations  thought  to  be  important  and  then  solve 
this  "truncated"  CI  problem.   If  so,  the  full  CI  features  are  no  longer 
generally  true.   In  particular,  the  explicit  combinations  [Z, .]   are  impor- 
tant to  the  results  obtained, 

A  second  problem  encountered,  the  choice  of  basis  fx  1 ,  is  a  universal 
one  for  all  methods  v/ith  the  exception  of  purely  numerical  solutions.   In 
some  ways,  however,  the  selection  is  more  crucial  for  CI  studies  than  in 
others,  such  as  the  HF.   The  reason  for  this  is  two-fold.   First,  contrary 
to  a   HF  calculation,  if  the  CI  procedure  is  employed,  one  is  generally 
interested  in  excited  states  as  well  as  the  ground  state.   Tliis  usually 
requires  a  basis  that  can  do  both  adequately,  a  task  v/hich  is  exceedingly 
difficult  to  manage.   Also,  there  is  little  or  no  "chemical  intuition" 
which  can  suggest  orbital  exponents  that  will  allow  a  satisfactory  de- 
scription of  excited  states.   Of  course,  if  the  basis  set  v/ere  large  enough, 

then  the  choice  of  the  fx  ]  v7ould  not  be  as  crucial,  and  this  is  a  fairly 

(J, 
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reasonable  alternative  in  HF  theory,  but  in  CI  one  is  again  led  to  the 
problem  with  the  number  of  possible  configurations. 

One  has  essentially  attempted  to  resolve  this  dilemma  with  two  meth- 
ods, both  used  separately  and  in  conjunction.   These  are  orbital  exponent 
optimization,  which  is  most  effective  when  undertaken  for  each  state  in- 
dependently, and  natural  orbital  iterations.   The  first  of  these  is  self- 
explanatory,  although  still  largely  prohibitive  due  to  the  extreme  cost. 
The  second  is  a  potent  tool  for  improving  less  than  full  CI  calculations 
by  exploiting  the  flexibility  in  the  selections  of  the  linear  combinations 
{5.}.   This  technique  makes  use  of  the  concept  of  "natural  spin  orbitals," 
introduced  by  Lowdin  [9],  vzhich  are  defined  as  the  eigenvectors  of  the 
first  order  reduced  density  matrix, 

Y(xjx')  =  [^"(x  x^.,.x  )  Y(x'x„...x  )dx.,...dx      (1.9) 
11    J     1^    n     1^    n/     n 

in  the  sense  of 

jY(xJxpcpj^(xpdx|   =  VkK^    '  ^^-^^^ 

where  n^  is  the  occupation  number  associated  v/ith  the  natural  spin  orbital 
cp,  .   From  a  theorem  due  to  Schmidt  and  later  Coleman  [lO],  it  is  found  that 
this  set  of  orbitals  has  a  certain  optimum  propertry.   That  is,  for  a  full 
CI  wavefunction,  the  selection  of  the  natural  orbitals  with  largest  occu- 
pation numbers  as  the  spin  orbitals  from  which  the  configurations  are 
constructed  results  in  a  much  shorter  approximate  expansion  which  has  the 
property  of  being  in  maximum  possible  coincidence  with  the  full  CI  [11], 
Natural  orbital  iterations  further  utilize  this  feature,  except  now  per- 
taining to  a  truncated  CI.   Since  such  an  approximation  is  not  invariant 
to  the  choice  of  the  {§.},  by  obtaining  a  set  of  approximate  natural 


orbitals  from  the  density  matrix  associated  with  this  CI  wavefuuction,  one 
can  select  a  set  via  an  occupation  number  criterion  from  v;hich  another 
trial  CI  wavef unction  is  generated.   This  set  is  expected  to  be  an  improve- 
ment over  the  previous  {5.}.   Consequently,  by  continuing  the  process  and 
assuming  convergence,  the  mixing  of  the  basis  functions  {x  }  is  thusly 
"optimized"  [12,13].   This  technique  is  very  useful,  but  it  also  has  de- 
ficiencies, one  of  v;hich  is  the  dependence  on  the  initial  selection  of 
configurations . 

Having  outlined  the  essential  components  of  the  two  most  important 
quantum  chemical  methods,  given  a  satisfactory  basis  set,  the  predominant 
problem  that  emerges  in  progressing  beyond  the  Hartree-Fock  approximation, 
as  well  as  utilizing  configuration  interaction  to  maximum  efficiency,  lies 
in  the  exorbitant  numbers  of  configurations  that  need  to  be  considered. 
The  treatment  of  this  standard  problem  in  various  contexts  is  the  objective 
of  this  investigation.   By  appealing  to  the  partitioning  technique  [14,15] 
and  perturbation  theory,  we  shall  develop  an  approach  to  this  problem  that 
differs  from  the  natural  orbital  methods,  although  it  is  entirely  com- 
plementary.  This  process,  the  "reduced  partitioning  procedure,"  shows 
promise  of  being  useful  in  a  variety  of  applications  within  the  configu- 
ration interaction  domain. 

The  basic  idea  of  the  procedure  is  to  reduce  the  dimension  of  a  CI 
problem  by  grouping  N  configurations  into  M  functions,  such  that  the  energy 
can  be  obtained  to  a  desired  numerical  accuracy  even  though  M  «  N.   The 
main  distinction  between  this  "reduction"  procedure  and  the  natural  orbi- 
tal methods  is  that  the  M  functions  obtained  arc  linear  combinations  of 
the  previous  coi^f igurations ,  whereas  a  natural  orbital  "reduction"  replaces 
one  set  of  configurations  with  a  second.   Thus,  since  each  technique  em- 
phasizes different  aspects  of  the  problem,  they  can  be  used  in  conjunction 
for  additional  flexibility. 
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The  method  employed  to  achieve  this  result  is  essentially  perturbation 
theory,  but  in  a  mora  general  framework  than  is  customary.   Traditionally, 
the  feasibility  of  applying  perturbation  theory  to  problems  has  depended 
on  three  elements:   the  existence  of  an  exact  solution  to  an  unperturbed 
Hamiltonian;  the  convergence  of  a  linear  perturbation  expansion;  and, 
due  to  inherent  computational  problems,  the  suitability  of  first  and  sec- 
ond order  solutions.   These  are  impractical  restrictions  for  arbitrary, 
many-electron  systems,  however,  and  preferably  variational  principles 
should  be  incorporated  with  perturbation  techniques  to  give  more  powerful 
methods  for  the  calculation  of  molecular  properties  and  energies.   The 
perturbation  approach  we  have  fonnulated  includes  the  possibility  of  using 
any  approximate  v/avefunction  as  the  unperturbed  solution,  guarantees  con- 
vergence via  non- linear  summations,  and  allows  higher  order  solutions  to 
be  easily  obtained. 

The   reduced  partitioning  procedure   (RPF)  will  have  several  advan- 
tages over  the  more  standard  methods.   For  instance,  a  full  CI  calcula- 
tion is  seldom  soluble  since  the  diagonalization  of  large  Hamiltonian 
matrices,  extremely  time  consuming  at  best,  is  usually  impossible  due  to 
near  linear  dependencies.   Although  the  reduction  process  would  generally 
require  that  all  the  matrix  elements  are  constructed,  it  would  oiily  require 
a  diagonalization  of  an  M  x  M  Hamiltonian  matrix  for  M  «  N,  thereby  pro- 
viding a  feasible  alternative  to  the  full  CI  solution.   Even  though  a  full 
CI  may  not  be  the  objective,  one  may  be  interested  in  certain  intermediate 
solutions  v/hich  would  give  a  substantial  improvement  over  a  Hartree-Fock 
or   truncated   CI  v;ave function.   In  this  respect,  also,  the  RPP  suggests 
some  possibilities.   One  of  these,  the  first  order  solution,  yields  a 
dramatic  energy  decrease  while  possessing  important  computational  simpli- 
fications.  This  first  order  solution  can  also  be  connected  with  the 
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"geometric"  approximation  and  the  method  of  "steepest  descent."   For 
excited  states,  the  RPP  uill  provide  two  alternative  approaches  which 
have  different  advantages.   One  of  these  allows  a  selective  treatment  of 
the  particular  state  of  interest.   The  question  of  properties  other  than 
the  energy  is  also  a  significant  one,  since  an  ordinary  perturbation 
treatment  of  second  order  properties  from  a  CI  perspective  would  require 
a  set  of  N-1  excited  states,  or,  equivalently,  the  inversion  of  an  N-1  x 
N-1  matrix.   Consequently,  here  also,  a  systematic  reduction  of  the 
space  is  v/arranted,  as  this  v/ould  permit  a  reasonable  approach  within  the 
CI  framework. 

After  presenting  background  material  on  the  partitioning  technique, 
perturbation  theory,  and  Fade  approximants  in  the  next  two  chapters,  the 
general  formulation  of  perturbation  theory  and  the  reduced  partitioning 
procedure  will  be  developed  in  Chapter  IV,  with  the  discussion  of  excited 
states  following  in  Chapter  V.   In  the  final  chapter,  numerical  applica- 
tions to  the  H„  molecule  and  the  HeH  molecular  ion  will  be  described. 


CHAPTER  II.   THE  PARTITIONING  TECHNIQUE  AND 
PERTURBATION  THIiORY 


The  partitioning  technique,  as  developed  by  LbVdin  [14,15],  has  been 
shown  to  be  of  primary  importance  in  elucidating  the  inter-relationship 
between  perturbation  theory  aiid  the  eigenvalue  equation  as  vrell  as  leading 
to  a  theory  of  upper  and  lower  bounds  to  energy  eigenvalues.   In  the  fol- 
lowing, the  partitioning  technique  will  be  briefly  presented  from  which 
the  resolvent  formulation  of  perturbation  theory  will  arise  as  a  logical 
consequence.   In  addition,  the  important  idea  of  inner  and  outer  projec- 
tion will  be  described  in  the  last  section. 

II. 1.  The  Partitioning  Technique 
Tlie  paramount  objective  of  molecular  quantum  mechanics  is  the  solution 
of  the  stationary  state  Schrodinger  equation, 

HY  =  EY  ,  (2.1) 

where  H  is  the  Hamiltonian  for  the  system  and  Y  an  eigenf unction.   If  we 

assume  an  expansion  for  Y  in  terms  of  an  orthonormal  n-particle  basis 

ff )  =   f^f„...f  ...),  then  instead  of  (2.1)  we  obtain  the  matrix  form  of 
'      '  i  Z    n 

the  eigenvalue  problem 

IIKE  =  (E  IE  (2.2) 

where  the  expansion  coefficients  for  each  state  are  grouped  into  d 
((I  C^  ...C   ...).   If  Jf  is  complete,  then  (2.1)  and  (2.2)  are  equivalent. 
For  a  finite  basis,  assumed  to  be  of  the  order  W-1 ,  the  eigenvalues 
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(E  E  . ,.E  )  are  respectively  upper  bounds  to  the  exact  solutions.   By 
o  i    n 

defining  M  =  H  -  £•!,  (2.2)  can  be  written 


na   =  0 


(2.3) 


for  Ml  =  <jE|  M|ff), 

If  we  noxi7  consider  |lf)to  be  partitioned  into  two  subsets  |A)nndj(b) 

of  order  n   and  n,  ,  respectively,  then  it  follows  that 
a      b 


and 


1^ 


^3    MI  , 

aa    ab 

ba    bb 


(E   = 


^nO   ^al 

So     Si 


an 


bn 


From  (2.3)  and  (2.4),  we  have  the  system  of  equations 


(2.4a) 


(2.4b) 


n       C  .  +  H  ^  C  ,  .  =  (D 
aa  ai     ab  bi     a 


(2.5a) 


^a,  d  .  +  n,,ac ,  . 

ba   ai     bb   bi 


(2.5b) 


for  each  i.   By  solving  (2.5b)  for  <£.  .,    one  obtains 


(E  ,  .  =  -Ml,"/  Ml,   t  .  , 
bi      bb   ba   ai 


(2.6) 


as  long  as  tl         exists.   The  substitution  of  (2.6)  into  (2.5a)  gives  the 
expression 


aa     ab   bb   ba   ai     a 


(2.7) 


or,  equivalently,    if  written    in   terms    of   the   ilamiltonian  matrix, 


[(II       -  £1      )  +  m  ,  (ea, ,   -  in,,)"  ia,   ](t    .  =  o  .      (2.8) 

aa  aa  ab         bb  bb  ba      ai  a 
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Consequently,  by  defining  the  modified  Hamiltonian  matrix 

"»aa  =^'aa  -^"^^ab^^Sb"  ^'bb^'^^Ha'  (^.9) 

we  have  the  altered  eigenvalue  problem 

H   (C  .  =  C  .  e.  ,  (2.10) 

aa  ai     ai  i 

which  yields  n  roots  for  £  =  S.  =  E.. 
a  1    L 

The  choice  of  partitioning  is  completely  arbitrary,  and,  in  the  course 
of  this  work,  various  choices  will  be  made.   For  the  present,  however,  let 
us  restrict  ourselves  to  the  case  when  n  =1,  and  call  this  single  func- 

3. 

tion   f     =  CD    ,    with  C      7^   0.      From    (2.10),    we  have 
00  o 

e     =   H        =   H       +    H    ,(£3  -    e,,)'\l|,     ,  (2.11a) 

o  00  00  ob         bb  bb  bo 


or,    with  E°  =   H        =    (cp    IhIcd   ), 
o  00  ^o'     ' ^o 

£q   =   E°+    (co^|h|  b)  (£1  i^j^-    <ro|H|  to»-\fb|H|cp^)  .  (2. lib) 

This  implicit  function  of  £,  (2.11),  is  the  "bracketing  function"  which  is 

crucial  to  the  theory  of  upper  and  lower  bounds  [15].   The  terminology 

comes  from  the  fact  that  for  any  £  >  E  ,  then  £  <  E   v/hich  allows  one  to 

o        00 

"bracket"  the  eigenvalue  E  . 

o 

The  connection  between  (2.11)  and  Brillouin-Wigner  [16]  perturbation 
theory  could  be  obtained  directly  by  separating  the  inverse  matrix  into 
its  diagonal  and  non-diagonal  parts  and  then  expanding  with  respect  to  the 
latter  by  means  of  the  relation 

(A-  IB)'  =  A"  -I-  A"  DiA"  +  A"^  lBA\""^lBA'^-f-  (2.12) 

The  relationship  with  perturbation  theory,  however,  is  more  conve- 
niently presented  in  terms  of  the  operator  formulation  of  the  partitioning 
technique. 
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II. 2.  The  Partitionli-ig  Technique  In  Operator  Form 
In  order  to  discuss  the  operator  formulation  it  is  necessary  to  in- 
troduce the  pt"ojector 

which  projects  the  reference  function  cp   from  the  space  of  interest,   e*  Is 
idempotent  and  self-adjoint,  i.e., 

(?  ^  (5,      ^^  ^  ,  >  (2.14) 

f 

where  9'  indicates  the  Hermitian  adjoint.   In  addition,  the  projector  Tor 

the  "orthogonal  complement"  to  cp  ,  the  space  previously  spanned  by  j  Ei),  is 

P  =  1  -  0',  (2.13) 

2         t 

with  P  =  P  and  P   -  P  as  before.   Also,  since 

OP  =  p©'^  0,  (2.16) 

the  projectors  O*  and  P  are  said  to  be  "mutually  exclusive." 

In  order  to  treat  (2.1)  we  shall  also  need  the  definition 

H  =--  PHP  (2.17) 

V7hich  is  an  "outer  projection"  of  H  v/ith  respect  to  the  subspace  of  P 
(see  II. 5).   The  outer  projected  Hamiltonian  (2.17)  satisfies  the  relation 

vzhere  Y,  and  E   are  its  eigenfunctions  and  eigenvalues,  respectively. 

The  reference  function  cp  is  also  a  trivial  eigenfunction  of  H  with 

^o 

eigenvalue  zero.   All  the  other  eigenfunctions  {Y,  }  arc  restricted  to  be 
in  the  subspace  of  P,   From  a  simple  theorem  about  outer  projections,  it 
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is  easy  to  see  that  E   s  E   in  order,  vjhere  {E  )  arc  the  eigenvalues  of 
H  for  the  eigenf unctions  Y,  [15(PTX)], 

Let  us  now  introduce  the  complex  variable,  g.  ,  and  subtract  £  from 
both  sides  of  (2.18)  to  give 


(s  -  H)Yj^  =  (e  -  E^)^^  .  (2.iy) 

Using  the  fact  that  PY,  =  Y,  ,  and  defining  the  "reduced  resolvent," 

T  =  (S  -  U)"-^?  ,      •  (2.20) 


it  follows  that  (2.19)  becomes 

TYj^  =  (&-  \)"\  ,  (2.21) 

demonstrating  that  T  and  H  have  the  same  eigenfunctions .   Therefore,  the 
eigenvalue  problem  (2.21)  is  completely  equivalent  to  (2.18),  but  T  also 
contains  the  complex  variable  £  vihich  adds  to  the  flexibility  of  the 
problem.   It  can  be  easily  shovrn  that 

T  =  TP  •-=  PT  =  PTP  ,  (2.22) 

and  from  (2.20),  one  obtains  the  useful  relationship 

P(e  -  H)T  =  P  ,  (2.23) 

V7hich  constitutes  a  connecting  link  with  the  eigenvalue  problem  (2.1). 
Let  us  now  define  a  "trial"  wavefunction  Y^  which  is  dependent  on 
the  variable  S  by 

\  =  ^o-'-  V^o  •  (2.24) 

Since  T  Hep   is  in  the  orthogonal  complement  to  cp  ,  Yp  is  composed  of  the 
reference  function  and  a  term  from  the  complementary  space.   With  the 
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definition  (2.24)  we  also  have  the  "intermediate"  normalization  condition 
satisfied,  that  is 

The  intermediate  normalization  is  a  convenient  choice  to  make  since  it  is 
applicable  to  both  the  discrete  and  continuous  parts  of  the  spectrum.   As 
long  as  the  normalization  integral  (Yp|Y„)  exists,  Y   is  said  to  belong 
to  the  discrete  part.   Using  (2.23)  it  is  seen  that  the  trial  wavef unction 
satisfies  the  algebraic  identity 

P(H  -  £)Yg  =  0  (2.26) 

for  all  the  values  of  S,   Also,  by  using  the  fact  that  G- +  P  =  1 ,  we  have 

(H  -  £)Yg  =  cp^(£^  -  £)  (2.27) 

where  S   is  defined  as 
o 

^o  =  <9o|h|V  =  <cpjH+  HTH[cp^)  .  (2.28) 

Expression  (2.28)  is  the  operator  form  of  the  "bracketing  function"  that 

is  equivalent  to  (2.11).   Hence  Y   satisfies  an  inhomogeneous  equation 

which  becomes  identical  to  the  eigenfunction  in  (2.1)  v;hen  £  =  £  =  E  . 

o        o 

The  trial  function  (2.24)  and  the  expression  (2.28)  are  essential  to  the 
following  treatment  of  Brillouin-Wigner  and  Rayleigh-Schrodinger  pertur- 
bation theory. 

II. 3.  Brillouin-Wigner  Perturbation  Theory  [16] 
Assume  that  the  Ilamiltonian  for  the  system  can  be  separated  such  that 

H  =  H  +  V,  (2.29) 

where  V   is    a   perturbation  v/hich  may  or   may  not   be   small.      Then    if    the 
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previously  defined  reference  function  cp   is  assumed  to  be  the  unperturbed 
solution  V7ith  Y  the  eigenfunction  to  the  entire  Hamiltoiiian,  one  may 
arbitrarily  define 

Y  =  Wcp  '  (2.30) 

vjhere  W,  often  termed  the  "v/ave  operator,"  takes  cp   into  Y.   Agaiii,  as 
throughout  this  paper,  we  shall  require  the  intermediate  normalization 

<cp^(Y>  =  1.  (2.31) 

By  multiplying  the  eigenvalue  equation  (2.1)  by  cp   (the  complex  conjugate 
of  cp  )  and  integrating  we  have 

E   =  <cpjH|Y)  =  <cpjH^+v|Y>  =  E°+  <cpjvw|cp^>,        (2.32) 

giving  the  energy  difference  between  the  unperturbed  result  and  the  exact 

solution  as  an  expectation  value  over  cp  of  the  quantity  W.   As  yet  v/e 
still  have  no  realization  for  W;  but,  by  appealing  to  the  previous  defini- 
tions (2,24)  and  (2.28)  and  using  the  relationship  H  cp  =  E  cp  ,  we  obtain 

o^o    o^o 

Yg  =  (1  +  T^V)cp^  (2.33a) 

and 

^n  =  1^+  <^PnI^+  VT  V|^PJ-  (2.33b) 


For  the  case  of  £  =  S  =  E  ,  we  have  complete  agreement  with  (2.32)  for 

o   .  o  I       ^ 

W  =  1  +  T  V,  which  defines  the  primary  formulas  of  perturbation  theory. 

In  all  perturbation  treatments,  we  want  to  have  expressions  dependent 

on  H  and  cp  .   Thus  it  is  convenient  to  define  the  reduced  resolvent  T 
o     ^o  o 

associated  with  H   in  an  analogous  manner  to  (2.20),  i.e., 

T   =  (£  -  PH  ?)"■*"?.  (2.34) 
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T   can  be  related  to  T  by  means  of  the  closed  form  of  the  Identity  (2.12), 
o 

(A  •-  B)"-^  =  A""^  +  a"^B(A  -  B)"""",  (2.35) 


now  pertaining  to  operators  rather  thaii  matrices.   With  the  separation 
(2.29)  and  using  (2.35),  we  find  that  T  and  T   are  related  by 

T  =  T  +  T  VT,  (2.36) 

o    o 

allowing  (2.32)  to  be  written 

£  =  E°  +  (cp  Iv  +  VT  V  +  VT  VtU  ).  (2.37) 

o    o    ^o'       o      o   '  ^o 

By  repeated  use  of  (2.36),  an  infinite  geometric  expansion  for  T-,  in  the 
quantity  T  V  may  be  obtained, 

CO 

T^  =  S   (T  V)\  ,  (2.38) 

^   k=0   °    ° 

leading  to 

Y.  =  cp  +  E   (T  V)\      ,  (2.39a) 

c-   ^o   ,  ,   o   ^o 
k=l 

and 

£  =  E°  -1-  <cD  |V[E  (T  V)"   ]|cp  >.  (2.39b) 

O       O        O  '      1   r>  O  '   O 

k=0 

To  further  illustrate  the  connection  with  the  more  standard  inhomogeneous 
equation  approach  to  Brillouin-Wigner  (BW)  perturbation  theory,  we  can 
define 

^k  =  ^V)\  ^^•''°''^ 

and 

Vl  =  <^ol^l^Pk>  '  ^^'"^^^^ 
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then  once  S  =  £  =  E  ,  (2.39)  becomes 


k=l 
From  the  recursion  formula, 

Vl  =  ^V>\  '  (2.42) 

and  (2.23),  V7e  have  the  set  of  inhomogeneous  equations 

Consequently,  as  long  as  it  is  possible  to  solve  the  equatioiis  (2,43),  one 
can  construct  all  the  {en,  }  from  the  unperturbed  solution  cp  and  the  succes- 
sive evaluation  of  the  fi,  .,  's. 

k+1 

II. 4.  Rayleigh-Schrodinger  Perturbation  Theory  [17] 
In  most  applications  of  perturbation  theory,  the  Rayleigh-Schrodinger 
(RS)  type  has  been  used  since  it  generally  has  better  convergence  proper- 
ties than  the  BW  variety,  as  well  as  not  yielding  equations  that  are  im- 
plicit functions  of  the  energy,  S.   This  latter  simplification  is  a  conse- 
quence of  the  unperturbed  energy  appearing  in  the  RS  reduced  resolvent 
rather  than  the  S  as  occurs  in  T  .   Instead  of  expanding  the  fundamental 
resolvent  T„  in  terms  of  V  as  before,  we  can  instead  consider 

V  =  V  -  (e  -  E*^)  (2.44) 

o 

as  the  quantity  to  be  used  in  the  expansion.   If  so,  from  (2.20)  and 
(2.29)  we  find  that 

T  =  R  +  R  V'T  (2.45) 

o    o  ' 

with  the  RS  reduced  resolvent  defined  as 
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R  =  (E°  -  ril  P)"^P.  (2.45) 

o     o     o 

In  analogy  v^^il;h  the  BVJ  case,  the  resolvent  and  wave  operator  become 

k 
T  =  E   (R  V)  R   ,  (2.47a) 

k=0  ° 
and 

k 
W  =  1  +  E   (R  V)  R  V.  (2.47b) 

k=0  ° 
The  connection  with  the  conventional  inhomogeneous  equations  of  RS 
theory  is  somewhat  more  complicated  than  in  the  BW  case  due  to  the  remain- 
der terras  which  arise  in  each  order  from  the  S  =  E  in  (2,44).   However, 
from  the  definitions  (2.30)  and  (2.32)  with  (2.47),  it  follows  that 

CO 

Y^  =  §  +  E   (R  V')^$  (2.48a) 

E    o   ,  1   o     o 

k=l 

and 

oo 

E  =  E°  +  E   <$  !V(R  V')^^U  >.  (2.48b) 

o    ,  -   o'    o     '  o 
k=0 

If  the  expressions  (2.48)  arc  arranged  after  powers  of  V  in  the  normal 

perturbation  series  expansions, 

00 

Y^  =  §  +  E  $,  (2.49a) 

E    o   ,  -,  k 
k=l 

and 

CO 

E  -  E*^  +  E   G,  ,  (2.49b) 

°   k=l  ^ 

then  the  lower  order  terms  may  be  identified  as 


$,  =  R  V§  (2.50) 

loo 

$_  =  R  (V  -  6, )R  va  =  R  (V  -  f..)^l. 
2  o      loo    o      11 
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and 


^1    =    <^ol^l^o> 

€„   =    <§    IvU,  )   =  <$    IVR  VU    >                                                        (2.51) 

2            o''l  o'o'o 

^4  =    <^JvhJ3>   =^  <^|VR,(V-€,)R^(V-6PR^V|$^> 


•€„<$    JVR^U    > 
2     o'      o    '    o 


Eventually   the   pattern  emerges   thtit   the   RS   perturbation  wavefunctions   $, 
can  be  defined  recursively   as 

k 

$,     ,    =  R    (V   -   €,)5      -   E     R  £      J,  (2.52) 

k+1  o  Ik        ._-,    o  k-^  X 

with  the  RS  energies  « imply 

^1  =  <^ol^l\>-  (2.53) 

These  are  the  usual  working  definitions  of  the  RS  quantities.   Using  the 
intermediate  normalization  condition,  (2.52)  is  easily  framed  into  the 
inhomogeneous  equation  form  by  means  of  the  simple  relation 

P(E°-  H  )R  =  R  rE°  -  H  )P  =  P  (2.54) 

000000 

that  is  obtained  analogously  to  (2.23). 

II. 5.  Inner  and  Outer  Projections 
Another  important  technique  which  is  often  used  in  connection  with 
perturbation  theory  and  variational  studies  is  the  method  of  "inner"  and 
"outer"  projections.   Since  we  shall  have  several  occasions  for  their  use, 
it  will  be  worthwhile  to  briefly  mention  some  of  their  properties. 
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An  outer  projection  of  an  operator  A  is  defined  as 

A  =  e-^e-^  (2,55) 

for  0     a  projector  characteristic  of  the  space  spanned  by  the  basis  [if). 
For  a  basis  of  either  finite  or  infinite  order,  6*^  may  be  written  explic- 
itly 

&^   =  |lf)<ff|iE)"\if|,  (2.56) 

2  t 

which  is  easily  seen  to  satisfy  the  relations  &     =  &     and  &     =  &  .      In 

the  event  that  the  infinite  basis  is  complete,  then  0'  is  simply  the 
identity  operator.   For  a  moment  let  us  define  A  as  the  Hamiltonian  oper- 
ator for  the  system  and  assume  that 

\   =  |ff)\  (2.57) 

is  an  eigenf unction  to  H.   Then  from  the  eigenvalue  relation 

"\=\\'  (2.58) 

we  have 

|if)«if|ff)'^^^(tj^  -  E^a^)   =  0.  (2.59) 

Since  it  is  assumed  that  the  basis  [if)  is  linearly  independent, 

v/here  the  metric  matrix  A^=  (^fju)  vjas  introduced.   Therefore,  the  outer 
projected  Hamiltonian  H  has  an  exactly  soluble  eigenvalue  equation,  which 
is  identical  to  (2.1)  if  |if)  is  complete  and  corresponds  to  the  varia- 
tional solution  if  the  basis  is  finite.   In  the  latter  case,  the  set  of 
eigenvalues  {E  }  are  guaranteed  to  be  upper  bounds  to  the  exact  results 

K. 
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in  order,  since  diagonalization  of  11   insures  that  each  Y,  is  orthogonal 
and  noninteracting  with  respect  to  those  of  lower  energy  [13,19],   A 
Hamiltonian  of  the  form  H  is  convenient  since  it  allows  one  to  formulate 
quantum  chemical  problems  in  a  consistent  fashion  even  though  an  approxi- 
mate solution  may  be  involved. 

If  F  and  G  are  two  Herraitian  operators  v/hich  fulfill  the  condition 
that 

<xIf|x)  >  <x|g|x>  (2.61) 

for   all   possible  X   in   the   coiranon  domain  of  F   and   G,    then   it   can  be   said 
that 

F  >  G      .  (2.62) 

Consequently,  if  F  >  0,  F  is  positive  definite. 

Now  let  us  require  that  the  operator  A  be  positive  definite.   If  so, 
we  can  define  the  inner  projection  of  A  as 

A'  =  A^6'^A^  (2,63) 

Since  0  s^  <3f^  ^  1,  A'  satisfies  the  operator  inequality 

0  <  A'  -s  A.  (2.64) 

By  making  the  transformation 

|lli  )  =  A^  [ff  )  (2.65a) 

and 

|s>=A"^|if),  (2.65b) 

we  obtain  tv;o  alternative  forms  for  A',  namely  the  Bazley  projection  [20], 


23 

A'    =    1  5i)  <ni|  A'^ha)'^  <Si|  ,  (2.66a) 


and   the  Aronszajn   fomi  [21], 

A-  =  Ajs)  <s|A|s>'\ffi|A.  (2.66b) 

The  inequality  (2.64)  can  only  be  concluded  for  A  >  0,  but  expressions 
(2.66)  will  still  converge  to  A  in  the  limit  of  a  complete  basis,  consti- 
tuting a  reasonable  approximation  eve:i  though  the  operator  of  interest 
has  no  unique  sign. 

One  important  example  of  the  use  of  an  inner  projection  is  found  in 
the  even  order  terms  of  a  perturbation  series  expansion.   For  instance, 

in  the  RS  case  for  R  <  0,  from  (2.46)  and  (2.51), 

o 

6.   =    <$    IVR  Vl§    ).  (2.67) 

2  o'       o    '    o 

With   P|[li)=    I  Hi),     A  =    (-R   ),    and    (2.64),    v/e  have   R'    >  R      giving 

^2   =    <*ol^''o''l^>   -   ^2'  ^^-^^^ 

which  gives  an  alternate  derivation  of  the  Hylleraas  variational  princi- 
ple [22,23].   It  can  be  similarly  shovm  that  as  long  as  $  ,$^,...,§   ^ 
are  exactly  knowi  [24], 

for  any  n,  a  result  obtained  by  Scherr  and  Knight  in  a  different  way  [25], 
(See  also  [26].) 


CHAPTER  III,  FADE  APPROXIMANTS 


In  all  applications  of  perturbation  theory  one  is  concerned  wii:h 
various  series  expansions  for  functions  or  energies  which  quite  often 
converge  slowly  or  not  at  all.   In  particular,  in  modern  perturbation 
theory  the  restriction  that  the  perturbation  be  small  has  been  esseuLially 
eliminated,  often  leading  to  very  serious  convergence  difficulties  that 
require  more  powerful  summation  techniques  in  order  to  extract  useful 
information.   One  possible  approach  to  this  type  of  problem  that  has  h^en 
sho^'m  to  be  fruitful  is  the  Pade  approximant  summation  procedure  [27], 

Although  all  the  convergence  properties  of  Pade  approximants  are  not 
yet  knotim  [23,29],  in  the  cases  vjhere  mathematical  justification  is  avail- 
able, the  range  of  convergence  for  a  Pade  approximant  has  been  found  to  be 
vastly  superior  to  that  of  an  ordinary  power  series.   In  numerous  other 
practical  applications  where  a  rigorous  mathematical  proof  is  lacking,  the 
results  seem  to  indicate  that  the  range  of  convergence  is  actually  greatly 
in  excess  of  that  for  the  situations  in  which  it  has  been  proved.   In  the 
following,  the  basic  concept  and  definitions  will  be  briefly  presented 
including  a  convenient  x'elationship  with  the  inner  projection  technique 
which  will  subsequently  bo  of  use. 

III.l.  Derivation  and  Definitions 
Consider  a  function  f(z)  which  has  a  formal  series  expansion 

00 

f(z)  =  E  a.z/.  (3.1) 

The  [N,M]  Pade  approximant  to  this  series  is  defined  as  the  ratio  of  two 
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polynomials   P(z)  and  Q(;5) 

[N,M]  =  P(z)/Q(E),  (3.2) 

which  are  respectively  of  degree  M  and  N  and  coincide  with  the  scries 
expansion  (3.1)  through  the  (M  +  N)    power.   By  requiring  the  latter  prop- 
erty, the  coefficients  of  P(z)  and  Q(z)  are  determined  by  the  condition 
that 

f(z)Q(z)  -  P(z)  -^   0(z'"+  ^+  S  (3.3) 

with  0(z         )  indicating  zero  up  to  terms  of  order  greater  than  M  +  N, 
To  prove  that  the  ratio  (3.2)  gives  a  unique  [n,M]  approximant  assume  that 
another  ratio  (call  it  P'/Q')  is  possible,  then 

f(z)Q(z)  -  P(z)  =  0(z"'+  ""^  S,  (3.4a) 

f(z)Q'(z)  -  P'(z)  =  0(z"'^  ™"^  ^.  (3.4b) 

Left  multiplying  (3.4a)  by  Q'  and  (3.4b)  by  Q,  the  first  terms  may  be 
eliminated  to  give 

P(z)Q'(z)  -  Q  (z)P'(z)  =  0(z™"^  ""  ''■  S.  (3.5) 

Since  PQ'  and  QP '  cannot  be  of  higher  order  than  M  +  N,  to  that  order  we 
have 

P/Q  =  p7Q'  (3.6) 

which  establishes  the  uniqueness  of  the  Pad/  approximant  [30j. 
If  we  write 

^      k 
P(z)  =  E   b  z  (3.7a) 

k=0  ^ 
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and 

N      k 
Q(z)  =  S   c  z    ,  (3.7b) 

k=0  ^ 

then  be  performing  the  multiplication  in  (3,3)  and  equating  like  powers  of 
z,  one  gets  the  set  of  equations 


a  c     =  b 
o  o  o 


a  c     4-   a   c,    =  b, 
o   o  o   1  1 


ac     +a     ^c,   +    ...+   ac      =b 
n  o  n-1    1  o  n  n 

(3.8) 


ac+a      ,c.    +    ...+    a        c      =b 
m  o  ra-1    1  m-n  n  m 


a,c     +ac+...+a  -c     =0 

mt-l    o  ml  ra-nfl   n 


a       c+a  ,c,   +    ...+   ac     =0. 

mfn  o  ml-n-1    1  ran 


Since  the  last  N  homogeneous  equations  involve  the  N  +  1  unknowns,  {c.}, 
there  is  oxi   infinity  of  solutions  possible;  but  since  the  Fade  approximant 
is  only  determined  to  within  a  constant  factor,  by  arbitrarily  choosing 
c   =1,  the  remaining  [c.]  can  be  obtained.   With  these  c.'s  the  fb.]  are 

o  1-^  11-^ 

determined  from  the  first  N  +  1  equations  of  (3.5)  and  thus  Q(z)  and  P(z) 
are  specified.   The  explicit  Fade  solution  to  (3.3)  is 
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[N,M]   = 


in-rH-1 


m-n|-2 


m 


m  , 

E     a.      z- 
J-n 


in-n-i-2 


m- 11-1-3 


m 
E 


m-tvi-l 


m-n+2 


m 


n 


j=n-l 


a.        -  z- 


in-ivl-2 

m-rH-3 

mj-l 
n-1 


rrvl-1 


iiH-2 


E     a.z- 


ml-1 


r,M-2 


iTH-n 


(3.9) 


subject  to  the  convention  that  a,  =  0  when  j  <  0,  and  aiiy  sums  which  have 
the  inital  element  of  higher  index  than  the  final  element  are  excluded. 

Another  form  of  (3.6)  that  is  particularly  convenient  for  computation 
is  the  inner  projection  form  [15(PTX  ),31,32].   (See  also  III. 2.)   If  we 
define 


(a.a    . . .a     ) 

11+1     H-N-1 


(3.10) 


and 


A 


1   ^1  1-1-1 


. .  .a 


H-N 


-1^' 


then  for  M  =  N  -l-  j,  (3.8)  can  be  written  [33] 


(3.11) 


[N,m-j]  =  E  a/  -,   ^^-'^^^^.,,  -   ^^^,,  ]"'..  .,,,      (3.12) 

JC/~-\) 

ere  the  poles  are  given  by  the   solutions  of  A\  .  ,  -  zA\  .  .   -  0. 

'   T+1      1+2  ' 


wh ^ ..-.^-   ,  -_  .-..  ,    -  .  .  . 

The  expression  (3.12)  could  be  shown  to  be  equivalent  to  (3. 'J)  by 
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standard  detei-minant  manipulations,  but  the  equality  can  be  established 
more  efficiently  by  expanding  the  inverse  in  (3,12)  by  (2.12);  then 

j  .    .    2N-1  J.       , 

+  "higher  order  terms"  .  (3.13) 


If  we  write 


-1 
>1 


b 


J+1 
Jt-2 


b. 


(3.14) 


jfN 


where  each  fb .  is  a  row  matrix  of   N    elements,  then  from  the  properties 
of  the  inverse 


\^Jl     =  \ji- 


(3.15) 


With  the  notation  1  ,  meaning  a  column  matrix  with  the  k   element  unity 
and  the  others  zero,  we  can  establish  the  relations 


H-1   j+k 


k  ' 


1  s  k  <;  N 


(3.16a) 


and 


j+2  k     j+k+1  ' 


(3.16b) 


from  which 


.V  -1   .   ^ 

'H-2'^j-H->   »  jfl   =  ^  y'ri' 


(3.17) 


and 


-1  -1  0 

j+1  "^  J4-2  jfV  jM  & 


(3.18) 
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From  (3,18),  (3,13)  finally  becomes 

2Nfj 
[N,Nl-j]  =  S  a  z  +  "higher  order  terms"  .  (3.19) 

X=0  ^ 

By  the  uniqueness  theorem,  it  necessarily  follov7s  that  (3,12)  is  equiva- 
lent to  (3,9)  [33]. 


III. 2.  Connection  with  Moment  Expansions 
and  Inner  Projections  [33] 


Let  us  consider  tV7o  operators  R  and  A  related  by 


R(z)  =  (1  -  zA(z)  )  '^ ,  (3,20) 


and  the  expectation  value  of  R  with  respect  to  a  normalized  function  ;|r, 

f(z)  =  <^lrlR(z)|t>.  (3.21) 

If  R(z)    is    expanded   by    (2,35),    it    follovjs    that 

f(z)    =  Z     z-^<^|a'^|0.  (3.22) 

Defining  a   =  {^\k    Ijf),  a   is  the  l       moment  of  A  and  (3.22)  is  the  moment 
expansion  of  f(z).   For  any  basis  |h),  an  inner  projected  approximation 
to  f (z)  can  be  given  in  the  Bazley  [20]  form  by 

f(z)  ~  (((J  h)  (Ii|  1  -  zA|  h)'^  (ra|0  .  (3,23) 

Since  the  basis  [  Qi)  is  arbitrary,  we  are  at  liberty  to  choose  the  functions 

I  Ih)  =  |i;fAiirA^,}f,..)  ,  (3.24) 

from  which  it  is  found  that 

f(z)  ~  aV   -  zA,  ]'^a  .  (3.25) 

o  o      i     o 

The  result  (3.25)  can  be  identified  from  (3.12)  as  the  [N,N-1]  Pads'' 
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approximant ,   Furthermore,  since  the  expansion  of  (3.20)  can  also  be 
written  as 

R(z)  =  1  +  zA(l  -  zA)'  ,  (3.26) 

it  is  possible  to  give  another  approximation 


f 


(z)  ~  a^  +  zi^'[A\^  -  zA\2  ]"  a^  ,  (3.27) 


which  is  seen  to  be  the  [NjN]  approximant.   Similarly,  the  higher  Fade' 
approximants , 

f(z)  ~  [N,Nf  j]  =  E   a^z^  +  z^\i^^(A^_^^  -  zA\^2)'^*^l'  ^^-^S) 

are  derived.   In  the  event  that  R  is  positive  definite,  from  the  inner 
projection  property  (2.64)  it  can  be  further  concluded  that 

[n,N]  s  f(z)  5  [n,N-1].    [33]  (3.29) 

In  this  case  the  moment  expansion  is  a  series  of  Stieltjes  and  many  addi- 
tional conclusions  may  be  reached.   For  other  operators  that  are  not  of 
a  definite  sign,  one  can  usually  only  treat  (3.28)  as  an  approximation 
to  f(z).   However,  since  in  the  subsequent  discussion,  the  [N,N-1]  Fade' 
approximant  will  be  identified  as  a  partitioned  eigenvalue  problem,  in 
this  instance,  it  necessarily  gives  an  upper  bound  for  any  finite  basis. 


CHAPTER  IV.  THE  REDUCED  PARTITIONING  PROCEDURE 

By  a  consideration  of  perturbation  theory  within  the  configuration 
interaction  (CI)  domain,  we  shall  develop  a  procedure  that  allows  an 
adequate  representation  for  a  CI  wavefunction  with  only  a  small  number  of 
perturbation  corrections  to  an  unperturbed  reference  function.   This  pro- 
cess, termed  the  reduced  partitioning  procedure,  leads  to  several  advan- 
tages over  the  more  standard  approaches,  one  of  which  is  found  in  the 
treatment  of  the  full  CI  problem.   The  theory  as  formulated  incorporates 
variational  principles  with  perturbation  methods  to  yield  a  scheme  that 
has  kno^im  convergence  properties  and  one  which  is  also  very  convenient 
from  the  computational  viewpoint.   In  connection  with  the  non- linear 
summation  techniques  invoked,  a  modification  of  ordinary  perturbation 
theory  is  obtained  from  which  an  alternate  set  of  perturbation  corrections 
to  the  energy  and  the  wavefunction  are  defined. 

IV. 1.  Non-Linear  Summation  of  Perturbation  Expansions 
In  order  to  approximately  solve  the  Schro'dinger  equation, 

HY  =  EY,  (4.1) 

via  perturbation  theory,  one  usually  assumes  a  separation  of  the  Hamilton- 
ian  operator 

H  =  H^  +  XV,  (4,2) 

with  H  chosen  to  have  an  eigenfunction  §   and  its  associated  eigenvalue 
o  o  " 
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6  ,  i.e. , 


o 


o  o    o  o  v+'j; 

Within  the  Rayleigh-Schrodinger  framework,  it  is  then  assumed  that  both 
Y  and  E  can  be  expanded  in  terms  of  an  "order"  parameter  X, 

2 

Y  =  $^  +\§^  +  \  $2  +  • • •  '  (4.4a) 

and 

E  =  £^  +  X6^  +  \t2  +  (4.4b) 

By  substituting  (4.4)  into  (4.1)  and  equating  the  various  "orders"  of  X, 
one  obtains  the  {5.}  as  solutions  of  the  inhomogeneous  equations 

(^o  -  V'n  =    (^  -  ^)Vl  -  io^n-/r  ^'-'^ 

with  the  perturbation  energies  [G  }  defined  as 

Recalling  the  definitions 


o    o  '  o  ' 


(4.7) 


and 


^o  =  ^^o  '  ^V^"^^'  (^•^> 

then  using  R  $  =0,  (2.54),  and  the  intermediate  normalization 

pI\)  =  \\)>  (4.9) 

(4.5)  may  be  conveniently  written  as 
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n-2 


f  =  R  (V  -  €,  )f  ,    -   T.     R  6   ,$,.  (4.10) 

n    o^     1'   n-1    ^^  o  n-^  Z 

If  the  expansions  (4.4)  are  to  yield  a  valid  solution,  they  must  con- 
verge to  the  exact  eigenfunction  and  eigenvalue  of  the  Hamiltonian,  H; 
but  this  may  not  be  fulfilled  for  many  choices  of  the  splitting  (4.2). 
Moreover,  it  is  seldom  possible  to  obtain  solutions  to  (4.10)  without  re- 
sorting to  some  approximate  method,  a  complication  which  tends  to  further 
obscure  the  validity  of  the  conventional  expansions.   To  increase  the 
flexibility  of  the  problem,  we  prefer  to  consider  the  more  general  finite 
expansions 

^(™)  =  $  +x$;  +  X.f:  +  ...  +  X""S'  ,  +  X"'^'       (4.11a) 
o    1     2  2  m-1       m 

and 

e^™^  =  £  +  \ei  +  \^a  +  ...  +  x™€'  +  x'^^e'  T  (4.11b) 

o     1      2  m        mfl 

for 

The  factors  [c.   },  which  may  be  dependent  on  the  truncation  order  m,  are 
at  our  disposal.   From  the  definitions  (4.6)  and  (4.12),  the  modified 
perturbation  energies  are  seen  to  be 

In  the  event  that  the  {c.   }  for  all  i  and  m  are  chosen  to  be  unity, 
standard  perturbation  theory  is  regained  leading  to  ^    =  Y  and  £    =  E 
at  infinite  order.   In  other  cases,  though,  we  have  an  additional  facility 
that  permits  us  to  insure  that  (4.11)  constitutes  a  valid  expansion  that 
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will  allow  f    and  £    to  approach  ijr  and  E  as  m  Increases  even  though  the 
ordinary  linear  perturbation  series  shows  no  convergence. 
The  Brillouin-Wigner  quantities  with 

T  =  (&  -  PH  P)"'^P,  (4.1.4) 

o         o 

are  defined  as 


^k+l  =  ^o^^)^^  (^•^5> 


and 


Vl  =  <^oI^I^Pk>'  (^•^^> 


where 


£  =  (cp^jV  +  VT(fi)v|cp^>.  (4.17) 

With  the  previous  modification,  we  similarly  obtain 

and 

^k+1  =  <^ol^I^>  =  4"^\+l  •  (^•^5> 

Besides  the  conventional  choice  of  c.  or  a.  as  unity  for  all  i,  we 

11 

have  the  option  of  determining  the  coefficients  from  the  variational  prin- 
ciple, or  more  generally  the  Lippmann-Schwinger  extremum  principle  [34], 
where  the  functions  {$ . }  or  {cp.}  are  chosen  as  the  basis  vectors.  The 
variational  principle  has  the  important  advantage  of  assuring  that  the 
modified  energy  expansion  (4.4b)  gives  a  converging  upper  bound  in  any 
order.  Another  possibility  is  to  determine  the  coefficients  by  a  Pade* 
summation  of  the  perturbation  energies  given  by  (4.6)  or  (4.16).   This 
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last  alternative  can  yield  numerous  different  solutions  of  varying  degrees 
of  desirability.   In  particular,  the  BW  energy  expansion,  ('^+.17),  for  a 
positive  definite  perturbation  and  for  £  less  than  the  first  excited  un- 
perturbed state,  is  a  series  of  Stieltjes.   From  this  it  may  be  concluded 
that    determining  the  [a.   }  from  the  [N,N]  or  [N,N-1]  approximants 
results  in  the  resultant  £    for  m  ^  1  yielding  respectively  upper  and 
lower  bounds  to  E  [33].  (See  III. 2.)   This  last  feature  can  also  be  re- 
lated to  the  Lippmann-Schwinger  principle  [35], 

IV. 2.  A  Model  Hamiltonian  for  a  Finite  Basis 

As  it  is  our  intention  to  study  energies  and  other  molecular  proper- 

I   \* 

ties  as  well  as  is  possible  within  the  n-particle  basis  |ff)  ,  we  are 

necessarily  limited  to  results  that  could  be  obtained  from  the  full  CI 
solution.   As  has  been  discussed,  this  in  itself  is  a  worthy  objective 
since  essentially  all  of  molecular  quantum  mechanics  makes  use  of  some 
basis  set,  and, as  such,  the  whole  field  can  rigorously  do  no  better  than 
a  full  CI  within  the  selected  basis. 

Moreover,  in  the  vast  majority  of  practical  applications  of  perturba- 
tion theory,  especially  if  higher  order  corrections  than  cp^  are  desired, 
the  functions  defined  by  (4.10)  or  (4.15)  must  be  calculated  in  an  approxi- 
mate manner,  usually  by  introducing  a  basis.   Once  the  perturbation  quanti- 
ties are  so  approximated,  another  question  arises,  namely  will  the  proper- 
ties of  these  approximate  solutions  be  those  of  the  exact  ones,  and  further, 
will  the  approximate  ones  converge  to  the  exact  ones  as  the  basis  tends  to 
completeness. 


"■The  basis  |ff)=  |f  f^...f  )  in  general  contains  all  the  configura- 
tions of  the  correct  s^m-unetry  wnich  can  be  constructed  from  a  suitably 
chosen  Slater  orbital  or  some  other  type  of  basis  set. 
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In  order  to  suitably  dispense  with  these  inconveniences  from  the  on- 
set, we  shall  explicitly  introduce  the  orthonormal  basis  |if)by  defining 
a  model  Hamiltonian  as  the  outer  projection  of  the  correct  Hamiltonian  for 
the  system,  then 

K=  &^m^  =    \S)  (s\e\s)(S\  (4.20) 

for 

e^  =  lff)<ff|.        '  (4.21) 

From  the  splitting  (4.2),  we  also  have 


K     =  e.H  e'^  (4.22) 

o    f  o  f 


and 


V=  e^^ve-^.  (4.23) 


If  we  now  partition  the  basis  Iff)  into  a  reference  function  cp  and  n  ele- 

.  o 

ments  required  to  be  in  the  orthogonal  complement  to  cp 


^o' 


|ff)=  |cp^,h>  .  (4.24) 


With  this  partitioning,  a  representation  for  the  projector  of  the  orthogo- 
nal complement  to  cp  is 


'^o 


P  =  1  -  \^)(tp\    =  \h)(h\.  (4.25) 


^o  '  '^o 


A  useful  relationship  between  the  functions  [  Ih)  and  an  outer  projected 


operator  Q     =  &  QSr     is  that 


<h|n  |h)  =  (tijnjlh).  (4.26) 


""The  partitioning  need  not  be  into  one  function  and  the  remaining  n, 
but  may  also  involve  a  linear  transformation  among  the  elements,  IfE).   The 
subsequent  discussion  and  particularly  (4.26)  will  still  follow. 
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By  using  the  result  of  (II. 5),  it  immediately  follows  that  the  model 
Hamiltonian  (4.20)  is  exactly  soluble,  i.e., 

KYi  =  E.  Y.,  (4.27) 

for 

with  (E  the  matrix  of  eigenvectors  to  the  Hamiltonian  matrix  defined  by 

|ffi).   In  addition,  (4.22)  defines  an  K  which  always  has  a  spectrum  of 

exact  unperturbed  solutions  whether  or  not  H  itself  has  a  set  of  eigen- 

functions.   Tliis  last  feature  eliminates  the  usual  perturbation  theory 

tyranny  of  limiting  one  to  a  very  few  possible  forms  for  H  ,  since  K     can 

o         o 

be  used  for  any  form  of  unperturbed  Hamiltonian  which  appears  to  be  desir- 
able.  In  the  event  that  one  actually  has  an  H   that  yields  a  spectrum  of 
eigenfunctions , 

then  by  choosing  Iff)=  |e  Q.  .  .  .Q    ),    it  follows  from  the  definitions  that 
'      '  o  1    n 

\\   =  ^\  (4.30) 

also.   Thus,  the  quantity  JC  can  be  viewed  as  another  generalization  of 
the  ordinary  theory. 

Relative  to  the  finite  basis  set  and  the  K  of  (4.22),  the  previous 
expressions  for  the  resolvents  R  and  T  become  exactly  [24] 

R  =|Bi)<h|£  -  K  I  ti)"  <Ih|  =  I  h>R  <lh|  (4.31a) 

o'      'oo'       ''o 

T^  =  |lh)(lh|e  -  3{:^|lh>"\[h|  =  I;ii)T  <Ih|,  (4.31b) 

where  the  relationships  p|lli)=  [h)  and  (4.26)  were  used.   If  the  functions 
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I  Ih)  are  eigenfunctions  of  K   ,    then  the  inverse  matrix  is  diagonal  giving 
the  more  customary  forms, 

R  =  2  -JS ^  (4.32a) 

°   k=l   £  -£° 
---  -,  o  k 

and 

T  =  S  -^ —  .  (4.32b) 

°   k=l  &-e^ 

Then,  subject  to  the  intermediate  normalization,  the  RS  expressions  are 

$   =  I  Ih)^   ,  (4.33a) 

n   '     n 

n-2 
«,^=  IR^[<nn|v  -6jlh>«^^.i  -  E  S^.^^^l  (4.33b> 


i=l 


and 


€^,  =  <$^|v|lh)*  ;  (4.33c) 

nf  1     o'  '    n 

and  for  the  BW, 

CD  =  |h)p  ,  (4.34a) 

^n    '    n 

le  =  T  (ni|vllh>l0   .  (4.34b) 

n    o   '  '     n-i 

and 

e  ,,  =  <cp  |v|  h)r3   .  (4.34c) 

nf  i     o '  '     n 

These  expressions  (4.33)  and  (4.34)  which  follow  from  the  model  Hamilton- 
ian  (4.20)  are  identical,  of  course,  to  those  obtained  from  the  Hylleraas, 
Scherr,  and  Knight  (HSK)  method  which  yields  approximate  perturbation  quan- 
tities for  the  exact  Hamiltonian  H  =  H  +  V  [22,24,25].   That  is,  instead 
of  solving  the  problem  (4.1)  approximately,  we  are  defining  a  Hamiltonian 
K  =  K    +y   whose  eigenfunctions  and  eigenvalues  can  be  exactly  determined. 
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Our  results,  however,  are  still  slightly  more  general,  since  our  JC   con- 

o 

tains  the  HSK  H   as  a  special  case, 
o 

From  (4.33)  and  (4.34)  the  elements  occurring  in  perturbation  theory- 
can  be  precisely  calculated,  and  the  coefficients  of  the  general  expan- 
sions (4.11)  and  (4.18)  can  be  chosen  to  guarantee  convergence  to  the  ax- 
act  result  for  JC,  ensuring  that  the  orthogonality  properties  of  the  eigen- 
function  to  an  Hermitian  operator  are  obeyed.   Furthermore,  as  [ff)  be- 
comes complete,  the  solutions  of  the  model  Hamiltonian  must  approach  the 
correct  eigenfunctions  and  eigenvalues  of  the  problem,  and  by  the  varia- 
tional theorem  always  from  above.   Similarly,  by  a  selective  choice  of  the 
coefficients  in  (4.11)  and  (4.18),  the  energy  to  any  order  of  solution  can 
also  be  made  an  upper  bound  --  and  sometimes  a  lower  bound  --  which,  con- 
trary to  much  of  modern  perturbation  theory,  allows  one  to  always  know  the 
relation  between  the  computed  result  and  experiment. 

IV. 3.  The  Fundamental  Equations  for  the 
Reduction  Process 

In  the  previous  section  we  defined  a  perfectly  general  K  as  a  multi- 
dimensional outer  projection.   We  shall  have  occasion  to  use  a  form  simi- 
lar to  this  in  the  discussion  of  excited  states,  but  for  the  present, 
another  form  arising  from  an  additional  outer  projection  relative  to  cp  , 

^o  =  ^^o^'       ®'=  l^o><^ol'  ^"^-^^^ 

* 
is  found  to  be  an  especially  convenient  choice.    Obviously,  this  ?/  sat- 
isfies (4.3)  as  well  as  allowing  us  to  introduce  any  approximate  wave- 
function  as  the  reference  function  in  a  perturbation  formulation.   From  the 
partitioning  (4.24),  the  definition  (4.22),  and 


-This  particular  form  for  V  has  been  previously  discussed  by  Musher 
[36]  and  Epstein  and  Karl  [37], 
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P|li)    =    |lh>,  (4.36) 


it   follows   that 

^o   =    l^o^Kl^><^l"ol^^^^'^o><^ol    =    l^o^^o^^ol  ("^-^^^ 

for  S   the  expectation  value  of  the  arbitrary  reference  function  cp  .   With 

this  %[   ,    we  have  the  definition 
o 

V  =  K  -  ,V  .  (4.38) 

o  • 

The  disadvantage  of  this  particular  /!/     is  that  it  leads  to  divergent  ex- 

o 

pansions  [36].   Thus  in  order  to  obtain  meaningful  results,  its  use  neces- 
sitates a  non-linear  summation  procedure  to  obtain  the  arbitrary  coeffi- 
ciants  in  (4.12)  or  (4.18). 

Another  feature  of  the  U     of  (4.35)  is  that  it  only  has  a  single  non- 
vanishing  eigenvalue.   This  results  in  the  R  of  (4.8)  and  T  of  (4.14) 

o  o 

assuming   the  simple   forms 


R     =  P/£      ;  T     =  P/e.  (4.39) 

o  o  o 

Restricting  ourselves  for  the  present  to  the  simpler  BW  expressions, 
in  the  following,  we  shall  deduce  a  reduction  procedure  by  appealing  to 
the  properties  of  a  moment  series  expansion  and  the  use  of  Pad^  approxi- 
mants .   The  final  closed  expression  can  be  identified  as  an  inner  projec- 
tion and,  as  such,  guarantees  an  upper  bound  in  all  orders.   From  the  de- 
fining equations  of  BW  perturbation  theory,  (2.37)  and  (2.40),  we  have 

£  =  <cn  IV  +  VT  V  +  VT  VT  V  +  ...  U  >  (4.40) 

^o '      o      o  o        '  ■  o 

and 
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By  using  (4.38)  and  (4.39),  with  P,V   =  0,  we  obtain 


o 


^=  ^o+   l/e<9o|KPK|cp^>+   l/£2(cp^|K(P3^)2|^^)+    ...      (4.42) 


and 


cp     =  1/  k(P3<:)''|cp  >.  (4.43) 


'o 


Then,    with   the   definitions 


cp^  =    (PJ{:)^|cp^)  =-|/li)<Ih|H||h)^"\lh|H|cp^)  (4.44) 

'hh 


=    ilh)    Hljjr^lb    ,  k  ^  1, 


and 


\-l   =  <^ol^l^k>  =  '^^"^hh^"''  ('►•'^5) 

the  l/£  moment  expansion  (4.42)  can  be  written  as 

^  ==  ^o  +  ^U   ^o  +  ^U^   '^l  +  •••  '  ^^'^'^^ 

where  the  previously  defined  BW  energies  are  &  =  CT  ^/c^    Consequently, 
given  the  Hamiltonian  matrix  in  the  complementary  space,  HI   ,  to  con- 
struct the  quantities  (4.44)  or  (4.45),  one  only  needs  to  perform  various 
matrix  multiplications. 

In  principle  it  should  be  possible  to  solve  (4.46)  as  an  implicit 
function  of  £,  and  if  the  convergence  is  sufficiently  rapid,  obtain  a  good 
answer  for  the  eigenvalue.   However,  this  is  impossible  for  two  important 
reasons.   First,  it  is  expected  that  the  series  (4.46)  should  be  strongly 
divergent.   For  example,  for  the  series  to  converge  the  eigenvalues  [^.3 
of  Hj^j^  should  be  such  that  -1  <  Q./e.  ^  +1    for  each  Q.    [38].   These  condi- 
tions are  unlikely  to  be  satisfied  for  molecular  problems.   A  second  com- 
plication is  that  for  any  truncation  of  (4.46)  there  is  necessarily  a 
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dependence  on  the  energy  zero-point.   However,  since  Fade  approximants  are 
known  to  have  better  convergence  and  existence  properties  than  a  linear 
series,  by  using  such  a  summation  technique  to  determine  the  coefficients 
in  the  elements  of  the  modified  energy  expansion  given  by  (4.19),  it  may 
be  possible  to  improve  the  situation.   In  particular,  by  applying  the 
[M,M-1]  Fade'' approximant  for  M  <  N,  (4.46)  can  be  written  as 

£^™^  =  £  +  l/,(ra)  [M,M-1],  (4.47) 

o     o 


where 


[M,M-1]  =  d'^iE      -  l/^(m)  E  )"\    ,  (4.48) 

oo      t      i      o' 

and 

CfT=  [p.a.....a.,      ,]  ,       i  =  0,1,2,..  .,in       (4.49) 
1   '-  1  1+1    i.+m-l-' 

With  some  slight  rearrangement  we  have 

£^"^  =  e  +  (r"^(&^™^E   -  E,)-\  .  (4.51) 

ooolo 

This  is  the  fundamental  formula  for  the  "reduction"  process,  but  as 
yet  it  is  not  obvious  if  the  two  previous  objections  are  eliminated.   Con- 
sequently, let  us  now  consider  the  functions  |jp)  defined  in  (4.44).   If  we 
augment  the  first  m  of  these  functions  by  cp  ,  and  then  assume  that  |cp  jp) 
is  a  reasonable  basis  for  an  approximate  solution  of  (4.1),  by  using  the 

fact  that  (cp  |cp  )  =  6    we  obtain  the  partitioned  form  for  the  secular 
^o ' ^k     ok 

equation  equivalent  to  (4.1)  as 

£^"^   =  £    +H      (fi^^^^A       -HI      )"^m       ,  (4.52) 

o  ocp     cp        cpcp        cpcp  cpo 

where  the  quantities  are  defined  as 


V  =<'?'"'!  ?'"*>■ 


H„-<f'"'|H|r<"'>, 
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(4.53a) 

(4.53b) 


and  £    is  the  eigenvalue  obtained  as  a  function  of  the  M  basis  functions 
||)^"^^).   From  (4.44)  and  (4.45)  it  also  follows  that 


H 


cpcp 


a 


'2    • 


o^ 


ra        nvt-1 


m 


2m- 1 


(4.54) 


and 


A 


cpcp 


..  a  . 

m-l 


a    T     a 

m-l        m 


'2m- 2 


(4.55) 


which  with  £    =  £    allows  the  identification  with  (4.51).   In  addi- 
cp 

tion,  the  intermediately  normalized  M   order  wavefunctions  are  given  by 


,  (^)  l~(ni)  s  .„  (m)  „     _  .-1 

t      =  9o+  la''  ')(£'  'e^  -  E^)     (j^. 


(4.56) 


Therefore,  the  summation  of  (4.46)  by  means  of  the  [M,M-1]  Fade*  approxi- 
mant  is  identical  to  solving  the  Hl-1  x  I'lfl  secular  equation  equivalent  to 
(4.52).  As  such,  besides  elucidating  the  nature  of  the  improved  result 
given  by  the  Fade  summation,  £    for  any  m  is  guaranteed  to  be  an  upper 
bound  to  the  full  solution;  and  the  result  obtained  from  a  variational 
calculation  must  be  independent  of  the  energy  zero-point.   In  the  event 
that  M  =  N,  we  have  simply  carried  out  a  linear  transformation  of  the 
original  problem  which  is  necessarily  equivalent  to  the  full  CI  solution. 
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The  expressions  (4.51)  and  (4.56)  are  also  seen  to  be  consistent  with 
the  definitions  (4.18)  and  (4.19)  v/here  the  arbitrary  coefficients  are 
obtained  via  the  variational  principle.   The  degree  of  success  of  the  re- 
duced partitioning  procedure,  however,  depends  on  the  [M,M-1]  Pade*  sxiin 
approximating  the  [n,N-1]  result  even  though  M  «  N,  since  then  £■ 
should  be  very  close  to  G    =  E,  and  hopefully,  to  a  desired  numerical 
accuracy.   If  so,  this  should  effectively  negate  the  necessity  of  solving 
enormous  secular  equations  that  would  arise  in  the  conventional  approach, 
as  well  as  leading  to  tractable  solutions  for  further  use. 

The  equivalence  of  the  [m,M-1]  Pad/  approximant  and  the  variational 
solution  in  the  basis  (cp  g)    )  could  have  been  concluded  in  a  slightly 
different  fashion  from  the  deduction  of  (III. 2).   This  result  shows  that 
a  moment  expansion  of  which  (4.46)  is  an  example  can  be  related  to  an 
inner  projection  which  in  turn  is  identified  as  an  [M,M-1]  Fade"' approxi- 
mant.  The  fact  that  it  is  an  upper  bound  follows  from  T  being  a  nega- 
tive definite  operator. 

Having  given  the  basic  equations  within  the  framework  of  BW  theory, 
it  may  be  asked  what  happens  in  the  RS  case.   Recalling  the  recursion 
expressions  (4.33)  and  (4.34),  and  by  using  (4.26),  (4.36),  (4.38),  and 
(4,39) ,  we  have 

\=  l'^>-k'    «'k=^/E  ^«hh"k-i  -^  A-r^^      (^-^7) 


X=l 


and 


%=    I^>l3k  '   "Pk  =  ^hh''\-l  '  (^-58) 

where  (4.58)  is  slightly  modified  from  (4.34)  to  be  consistent  with  (4.44), 
From  these  two  relations,  it  can  be  seen  that 
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av^  =  1/^  13^  (4.59) 


o 

o  o 

o  o 


This  can  be  written 

A  =  B  T  (4.60) 

for  A  =  {(K^(K„...(ii!    ),  IB  =  (I3t13o...|3  );  and,  where  it  is  also  ob- 
i   /     m  i   Z     m 

served  that  T  is  of  upper  triangular  form.  The  latter  fact  ensures  that 
for  any  order  m,  the  variational  solution  relative  to  either  basis  set  is 
the  same;  that  is, 

I  (™)       ■  I  5  (™)  \   (™)         I  (™)  \   (ni)  //  ^1  \ 

i    =  cp^  +  I  I    >«     =  cp^  +  h^    )*  (4.61) 

and 

£^™^  =  (^^"^IkU^"^.  (4.62) 

Although  the  variational  solutions  in  the  BW  and  RS  bases  are  equiv- 
alent, a  slightly  different  result  could  have  been  obtained   in  the  latter 
case  if  the  coefficients  {c .   }  were  determined  by  applying  a  [M,M-1] 
Fade"  approximant  to  the  series  (4.4b).   This  result  is  usually  only  slight- 
ly different  from  the  variational  solution  [24]. 

,IV.4.  The  Method  of  Steepest  Descent  and  the 

Geometric  Sumrule 

In  ordinary  perturbation  theory,  the  form  of  the  inhomogeneous  equa- 
tions for  the  perturbation  functions  is  dictated  by  the  requirement  that 
(4.4)  or  (4.40)  is  convergence.   In  the  case  of  a  divergence  expansion,  it 
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is  not  necessarily  obvious  why  the  set  of  perturbation  functions  should 
still  be  one  which  rapidly  exhausts  the  space  of  interest,  as  is  neces- 
sary for  the  reduction  process  to  be  effective.   To  attempt  some  justifi- 
cation  of  this  point,  we  appeal  to  the  method  of  "steepest  descent"   as 
employed  by  McWeeny  [40]  in  connection  with  SCF  studies. 

Given  an  arbitray  trial  function,  cp  =  cp  +  Sep,  where  it  is  assumed 
that  (cp  I  6cp)  =  0  and  &  =  (cp  |H|cp  ),  the  energy  increment  due  to  the 
first  order  variation  of  the  wavefunction  is 

(5cp|H   -   e  Jcp   )+    (cp   |h   -   e   |6cp)+   (6cp|H   -    £  J  6cp> 
5£  =  2__2 2 ° . 2 •      (4,63) 


<^olVn'^  +    (6cp|6cp> 


^o'  ^o 


If  [ h)    constitutes  a  basis  in  the  orthogonal  complement  to  qj  ,  we  can 
assume  an  expansion  for  6cp  of  the  form 

|5cp>  =    |ili><i^  (4.64) 

for  c  ^  the  expansion  coefficients.   If  we  neglect  the  second  order  term, 
(6cp|H  -  £  I  6cp)j  and  ask  for  the  direction  of  the  vector  in  the  linear 
space  corresponding  to  the  maximum  possible  change  in  ^S,,   we  obtain 

5(6£)  =  0  =  6a^  +  [(fh|H  -  £^|lh>  -  6e«^]  ,  (4.65) 

or,  that  a  is  proportional  to  (lli[H|cp  )  =  lb  as  defined  previously.  Con- 
sequently, the  basic  functional  form  for  the  first  order  correction  should 
be 

9^  =  llh>lb  .  (4.66) 

The  optimum  amount  of  cp^  to  mix  with  cp  ,  that  is  the  step  size  Tj,  can  be 


"Although  philosophically  similar,  the  terminology  "steepest  descent" 
usually  pertains  to  a  different  application  originally  discussed  by  Rie- 
mann  and  Debye  [39]. 
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determined  by  including  the  ignored  second  order  term  and  requiring  that 

for 

e(Tl)  =  (9^+  1>?jHlcp^  +  1^9^>/<cp^  +  1^^Jcp^+Tl9^).       (4.68) 

Then  we  have  the  result  that 

T\   =  (e(Tl)a^  -  o^y'^a^  '  (4.69) 

where  (4.66)  and  the  definitions  (4.45)  were  used.   This  result  gives  the 
first  order  equations  obtained  via  the  steepest  descent  procedure  as 

£^^^  =  fi  +  a  (e^^^a  -  o.y^a  (4.70) 

o    o      o    1    o  ' 

and 

*^^^  =  cp^  +  9i(&^^^a^  -  ^i^'^o^   .  (4.71) 

These  two  equations  can  be  identified  as  those  corresponding  to  (4.51)  and 
(4.56)  for  n  =  1,  or  the  [1,0]  Fade'' approximant  solutions. 

We  can  similarly  justify  the  higher  order  expressions  by  considering 
the  function  (4.71)  as  a  new  reference  function  and  repeating  the  proce- 
dure where  now  (cp  |  sO  =  0,  but  <ij;^   |  6\;i)  "f   0.   If  so,  by  again  excluding 
the  second  order  term,  we  have 

(^    u    > 

By  assuming  6il(  =  j  lli)  t  „  ,  we  find  that 
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where  the  constant  p,  indicates  a  mixing  of  some  amount  of  lb  ,   Thus,  the 
next  function  which  should  be  employed  to  augment  cp^  is  exactly  the  cp^  of 
(4.44).   By  a  continuation  of  this  argument,  the  remaining  cp's  can  be 
similarly  justified,  thereby  providing  an  alternative  explanation  from  a 
source  external  to  perturbation  theory  for  the  special  forms  of  (4.44). 

The  optimum  step  size  subject  to  cp  ,  cp^  ,  and  cp«  is  obtained  from  a  form 

(2)         —      - 

\|r    =  cp  +  Ucp,  +  T|pCppj  where  T|^  and  Tl„  are  variational ly  determined. 

This,  of  course,  is  the  [2,1]  Fade* approximant  of  (4.51)  and  (4.56). 
Another  approximation  could  have  been  obtained  by  assuming  the  \li  of 
(4.71)  to  be  a  new  reference  function  and  asking  for  the  steepest  change  in 
the  first  order  energy  relative  to  this  unperturbed  solution.  This  type 
of  procedure  has  been  employed  in  a  RS  framework  by  Dalgarno  and  Stewart 
[41]  and  Hirschfelder  [42,43].  For  the  BW  case,  compare  also  Goldhammer 
and  Feenberg  [44]  and  Young,  et  al.  [45]. 

Besides  being  directly  derivable  from  the  steepest  descent  procedure, 
the  first  order  solutions  (4.70)  and  (4.71)  have  some  other  properties 
which  make  them  particularly  significant.   For  example,  their  use  corre- 
sponds to  the  application  of  a  geometric  sum  rule  to  a  series,  in  this 
case  (4.46),  and  hence  is  often  referred  to  as  the  "geometric"  approxima- 
tion.  Despite  the  fact  that  such  an  approximation  is  seldom  warranted  by 
the  actual  terms  in  the  series,  it  has  been  used  with  a  great  deal  of  suc- 
cess in  a  variety  of  applications  [46,47,48,49].   This  approximation  -- 
recognized  as  the  [1,0]  approximant  --  has  been  previously  justified  on  a 
variational  basis  [50],  as  in  the  present  case,  where  it  rigorously  corre- 
sponds to  a  variational  solution  corrected  through  the  first  order  in  the 
wavefunction. 
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As  should  be  expected  from  the  connection  with  the  method  of  steep- 
est descent,  this  first  order  solution  gives  a  substantial  improvement 
over  the  unperturbed  result,  yielding,  in  the  caseswe  have  studied,  about 
sixty  to  seventy  per  cent  of  the  possible  energy  improvement  that  may  be 
gained  from  the  given  basis  set.   Furthermore,  the  geometric  solution 
also  has  an  important  simplifying  feature  when  the  reference  function  is 
chosen  to  be  the  two  most  common  types,  namely  the  Hartree-Fock-Roothaan 

(HFR)  or  the  CI  variety.   Let  us  consider  the  first  case  where  m      is 

^o 

chosen  to  be  an  HFR  wavefunction  and  assume  the  remainder  of  the  configu- 
rations are  generated  by  the  substitution  of  a  "virtual"  orbital.   If  we 
partition  |  h)  such  that  the  single,  double,  etc.  excitations  are  explic- 
itly introduced, 

|lh)  =  |lh^^)lh^2\..|h^">>,  (4.74) 

then  from  Brillouin's  theorem  it  follows  that 

^o  =  (%H^^^^)(^^^^N%)  (^-75) 

and 

o^  =  <cp^|H|lh^2^<tti(^V|lh^^^<{l/^)|H|cp^>.         (4.76) 

Consequently,  the  [1,0]  solutions  (4.70)  and  (4.71),  relative  to  an  IIFR 
function,  depend  solely  on  the  double  excitations.   The  CI  reference 
function  will  be  discussed  in  the  next  chapter  where  excited  states  are 
considered. 


CHAPTER  V.  EXCITED  STATES 

The  subject  of  excited  states  is  an  important  one  in  quantum  chem- 
istry, since  it  is  an  area  where  a  great  deal  of  information,  which  is 
not  readily  accessible  to  experiment,  can  be  obtained  from  theoretical 
consideration.   Compared  to  ground  states,  however,  the  difficulties 
associated  with  an  adequate  determination  of  excited  states  are  propor- 
tionately compounded.   In  particular,  it  is  expected  that  the  higher  ex- 
citations that  one  usually  neglects  in  a  CI  wavefunction  will  be  of  more 
significance  for  excited  states,  from  which  it  would  appear  that  a  reduc- 
tion process  should  be  of  even  more  importance.   In  the  following,  a 
super-matrix  formulation  for  excited  states,  which  is  a  generalization  of 
the  previous  theory,  will  be  presented.   From  a  consideration  of  these 
equations,  a  modification  is  made  which  allows  a  selective  treatment  of 
the  specific  state  of  interest.   This  last  result  will  also  suggest  a  con- 
ceptually appealing  approximate  treatment.   As  in  the  Hartree-Fock  case, 
the  first  order  solution  for  excited  states  is  found  to  have  a  computa- 
tionally convenient  property. 

V.l.  The  Super-Matrix  Formulation 
In  order  to  progress  to  excited  states,  it  is  necessary  to  consider 
a  multi-dimensional  reference  space   |^)  =  |x  X-i'-'X  )>  which  has  the 
characteristic  projector 


■''•"The  multi-dimensional  partitioning  technique  has  been  employed  in 
lower  bound  studies  by  Choi  [51]  and  Wilson  [52],  where  the  reference 
space  consisted  of  a  degenerate  set  of  functions. 
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\  =  lyXyl-  (5.1) 

Within  the  space  spanned  by  jff),  for 

Iff)  =  \-J(,^)  ,  '  (5.2) 


©*  is  the  identity  operator,  thus 

p  =  ©-^  -  e-^  =  |m)<ih|.  (5.3) 

Consider   an   eigenvalue   equation, 


30F.    =  Y.E.  (5.4) 

111  \    '    J 

for  some  state  Y.  with  eigenvalue  E..   Using  the  property  that 

K  =  &^  =  jce-^  =  ©^e-^  =  ©-^He*  ,  (5.5) 

(5.4)  can  be  rewritten  as 


3<:(0'  +  P)Y.  =  Y.E   .  (5.6) 

X      1    1  i 


From  (5.6)  the  two  equations, 


e-JKe-Y.     =©'Y.E.    -   e'KPY.  (5.7a) 

XXiXiiXi 


and 


PKOvY.    =  PE.Y.    -   PKPY.,  (5.7b) 

X   1       .       1    1  i'  ^  ^ 


are  obtained  by  left  multiplication  of  &     and  P,  respectively.   Solving 
(5.7b)  for  PY.,  it  follows  that 

PY.  =  Tg.J^e-^Y.  .  (5.8) 

where  T„  has  the  usual  definition 
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T^  =  (E.  -  PKP)'^?.  (5.9) 

E.     1. 

By  substituting  the  expression  (5.8)  for  PY. ,  (5.7a)  becomes 

e-^CK  +  KT  K)&  Y  =  e  Y.E  (j.lG) 

or  that  &   Y.  is  an  eigenfunction  to  the  operator  &   {K   +   JCT„  K)&      with 
XI  X        Ei    X 

eigenvalue  E..   (See  Lowdin  [15(PTIV)].)   Thus,  (5.10)  is  exactly  equiva- 
lent to  (5.4),  but  a  consideration  of  (5.10)  leads  us  to  a  perturbation 
treatment. 

With  the  definition  (5.1),  the  fact  that  Y.  =  |ff)C  .,  and  the  par- 
titioning (5.2),  (5.10)  in  matrix  form  is 

|y)<ylK+KTg.K|>?)<^Iff)C.  =  |^)(:3^1ff)C.E..     (5.11) 

For  l"^)  linearly  independent  and  <C .  =  (^^i),  this  requires  that 

hi 

This  gives  an  alternative  operator  derivation  of  (2.10),  which  was  pre- 
viously obtained  from  a  direct  partitioning  of  the  eigenvalue  equation. 
Obviously,  the  essential  difference  between  the  basic  equation  for  the 
multi-dimensional  case,  (5.12),  and  the  earlier  example  for  p  =  0  is  that 

now  we  cannot  remove  (£   to  get  a  direct  expression  for  the  energy  E.. 

Xi  1 

Instead,  we  must  solve  a  secular  equation  for  an  effective  Hamiltonian. 
This  consequence  ensures  that  the  resultant  multiple  state  solutions  will 
satisfy  the  requisite  orthogonality  and  noninteracting  properties.   Since 
the  expression  (5.12)  is  an  implicit  function  of  E.,  in  practice  we  must 
resort  to  the  more  general  form 


HI   (£)  (t  (&)  =  d  (e)e  ,  (5.13) 

XX     X       X 


53 

where  for  each  choice  of  £,  an  H   (£)  may  be  constructed,  from  which 

XX 

diagonalization  yields  p  +  1  roots,  {£  £^...£  }.   The  eigenvalues  [E . ) 
are  given  when  G  =  G,  =  E.  for  some  i. 
VJith  the  principal  definition, 

"^^^(-^^    =  (-^Ik+KT   K\^),  (5.14) 

we  can  proceed  as  previously  by  expanding  (5.14)  via  the  relation  T„   = 

T  +  T  KT,  to  obtain 
o    o 

H      (G)    =  H        +    <v|kt   K+  KT   KT   K+    ...|y).         (5.15) 
XX  XX  ^  '        o  o       o  '  *  ' 

With  the  definition 

€  ^(G)  =  <^|k(T^K)^+^|^),     k^O,  (5.16) 

(5.15)  can  be  written  as  an  expansion  in  terms  of  matrices 

e   (G)  =  ID   (G)  +  (£^(G)  +  €^(G)  +  €2(G)  +  (5.17) 

Now  let  us  assume  the  set  |^^)  are  the  solutions  to  an  unperturbed  pro- 
blem, that  is, 

for  t      a  diagonal  matrix  consisting  of  the  unperturbed  eigenvalues, 

[G  Gt..„G  ].   Then  the  expression  (5.17)  becomes  the  multi-dimensional 
o  1     p-^ 

equivalent  of  the  usual  BW  perturbation  expansions  which,  by  taking  to 
a  certain  order  and  then  solving  (5.12),  yields  the  perturbation  correc- 
tions for  each  of  the  p  +  1  roots.   This  type  of  treatment  for  degenerate 
RS  perturbation  theory  was  orginally  given  by  Van  Vleck  [53],   (See  also 
Hirschf elder  [54].)   For  a  general  JC   we  could  apply  various  treatments 
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to  (5.17)  to  obtain    series  of  non-linear  expressions  for  HI   ,  but  in 

XX 

the  interest  of  brevity,  we  shall  first  specify  a  particular  unperturbed 
Hamiltonian. 

Analogous  to  the  previous  single  partitioning  example,  we  shall 
select  for  the  unperturbed  Hamiltonian 


K 


X   X 


(5.19) 


which  is  easily  seen  to  satisfy  (5.18).   Using  the  spectral  expansion  and 
(5,3)  we  have,  just  as  before, 


T    =  p/e. 

o 


(5.20) 


Then  from  (5.16), 


k+1 


€^(r   =  i/.k+l  ^   =  l/^k+l  (>;Ik(PK)''+  |y>,        (5.21) 


where  the  [v  }  are  '    pendent  of  6.   By  applying  the  [m,M-1]  Fade'' approx- 
imant  to  (5.17)  we  obtain  to  order  m  the  super-matrix  expression  for  the 
reduced,  effective  HaT.iltonian  matrix. 


M^"''*   -  I    +  r"''(£ft    -  ft,  )'^r  . 

XX      o     o    o    1     o 


(5.22) 


The  super-matrices  [T ,  and  ft  are  defined  by 


(5.23a) 


and 


^,  =  [1-^r 


i^-i' 


.  r 


M-ra. 


'U 


Y, 


^1 


-1^ 


'i^-ra-l 


%1 


^Sl-m-1 


^. 


;^ra 


%2m-l 


X  =  0,1, 


(5.23b) 
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Inspection  of  (5.22)  indicates  that  the  super-matrix  product  results  in 
a  p+1  X  pfl  matrix.   A  somewhat  similar  kind  of  treatment  has  been  given 
in  [15(PTXIV)]. 

To  further  illustrate  the  connection  with  the  previous  theory,  we 
can  define  the  perturbation  corrections  relative  to  each  unperturbed 
state, 

cp^^^  =  (P3^)^|Xi>  =  |lh)H];;;\lh|H|x.)  =  |'h)lD;;j;^b.     (5.24) 

which  should  be  compared  to  (4.44).   Also,  we  have  the  analogous  sigma 
quantities 


.(i) 


(!)• 


t,„k 


a^-  =  (xj3C|cp-)=.Xh^ 


(5.25) 


and  for  i  5^  j , 


M>2)    _ 


,(j) 


.t,,k 


a^^'J^  =  <x.|K|cp^J^)  =  <x.|K|cp^^^)  =  -bJ-HlUj^lb.,      (5.26) 

which  are  generalizations  of  (4.45).   The  terms  a^     are  required  to 
describe  the  interactions  arising  from  the  different  states  |^).   With 
these  definitions,  the  elements  of  the  "^   matrices  are  seen  to  be 


\  -  ^\    "^k    •••  ^k   J 


(5.27) 


where  the  sigmas  are  grouped  into  column  matrices, 

.(0,i)' 


d 


(i) 


a, 


k 
(1,1) 


o, 


•(P,i) 

k 


(5.28) 


As  before,  we  can  identify  the  basic  reduction  equation  for  multiple 
states,  (5.22),  as  the  partitional  form  of  an  eigenvalue  equation,  but 
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this  time  relative  to  the  basis  [y^p    '£        ...g)   ),  for 

/'>  .  \.,['\['\..,,l''>).  (5.29) 

When  p  =  0  and  i  =  0,  the  fundamental  expression,  (5.22),  reduces  to  the 
single  dimension  case. 

V.2.  The  Selective  Reduction  Procedure 
The  advantage  of  a  reduction  process  based  on  (5.22)  is  that  all 
states  are  treated  equivalently,  giving  a  set  of  p  +  1  improved  states  at 
any  order  m.   This  is  also  a  disadvantage,  though,  since  in  order  to 
effect  this  result  one  must  use  (m  +  l)(p  +  1)  functions  in  the  computa- 
tion.  In  a  typical  CI  calculation,  even  though  we  expect  the  order  m  to 
be  low,  p  may  be  a  rather  large  number.   To  avoid  this  difficulty,  a  mod- 
ification of  this  method  will  be  presented,  which,  requiring  only  m  +  p  +  1 
functions,  permits  a  higher  order  treatment  for  one  specific  state. 

Although   it  is  possible  to  obtain  the  desired  result  directly  from 
a  consideration  of  the  super-matrix  equations,  for  illustrative  purposes 
it  is  convenient  to  approach  the  problem  from  a  different  perspective. 
Therefore,  let  us  consider  the  possibility  of  obtaining  a  wavefunction  for 
each  excited  state  from  a  perturbation  treatment  of  its  corresponding 
unperturbed  solution.   Then,  as  in  the  single  partitioning  example,  we 
assume  a  trial  function  of  the  form 

^e.  =  (^  +  ^e^^Xi,  (5.30) 

consisting  of  a  single  element  in  the  reference  space  |y),  and  a  terra 
from  the  orthogonal  complement,  [di).   Contrary  to  a  rigorous  single  parti- 
tioning, ilfg   cannot  become  the  correct  solution  Y.  when  &  =  S.  =  E.,  as 
long  as  Y.  contains  any  contributions  from  the  other  elements  of  |^). 
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However,  in  the  usual  application  of  the  theory,  the  set  [^  )  V7ill  be 
eigenvectors  of  a  relatively  small  CI  problem  and  will  be  necessarily 
mutually  orthogonal  and  noninteracting.   The  only  mixing  allowed  is 
found  in  the  Hamiltonian  matrix  elements  involving  the  remaining  (x . } 
and  the  terms  from  the  orthogonal  complement.   Consequently,  it  could  be 
expected  that  the  contribution  to  Y.  from  the  additional  (x . }  is  small. 
In  any  event,  we  can  obtain  the  approximate  expression 

E.  =  <x.|k|y^_>  =  (x^\K+KT^K\xr>,  (5.31) 

from  which  the  perturbation  functions  (5.24)  are  obtained,  and  the  m 
order  reduction  equations  analogous  to  (4.51)  and  (4.56)  are  derived; 

fiC^)  ^e?™^  .  gO,.  ^t(i)  (£E^^^-  Ei^h'\^^\  (5.32a) 

L         1         lOO         lo 

and 

^(m)^^(m)  ^^_^  |2)^'h(£E^'>  -  E^('))-'^^^\    (5.32b) 

One  way  to  retain  the  rigorous  upper  bound  properties  of  the  p  +  1 
approximate  solutions,   (5.32),  is  to  construct  the  Hamiltonian  and  over- 
lap matrices  relative  to  [ij;  \j[^...il(  }  and  to  solve  the  secular  problem. 
Then  one  regains  the  noninteracting  property  required  by  the  separation 
theorem  [18,19].   This  variation  is  not  very  flexible,  though,  since  the 
relationship  between  the  {x . }  and  {cp,   }  is  already  fixed.   A  superior 
alternative  can  be  deduced  by  proceeding  somewhat  differently.   Within  the 
space  of  interest,  and  subject  to  (5.19),  if  the  expression (5.14)  is  writ- 
ten in  more  detail,  we  have 


"The  solution  corresponding  to  x  >  the  ground  state,  is  still  ensured 
to  be  an  upper  bound  although  it  cannot  quite  converge  to  the  best  solution 
in  the  given  space  |ff)=  i^^   ,ni)  since  some  of  the  x-  s'^g  excluded. 
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HI     =  e  + 

XX    o 


^\h<^(^>^\o    "'oh^(^)^hl  •••  ^Ih^^^^^'hp 


^h^^^^^-ho 


>^(^>  '\o 


^h^^^^^'hp 


(5.33) 


for  Q(£)  =  (Si,,  -  HI   )   .   From  this  exact  expression,  (5.32)  is  seen 
to  correspond  to  a  reduction  process  applied  exclusively  to  the  diagonal 
elements.   Consequently,  assuming  the  perturbation  functions  so  obtained, 
[(p       cp        . ,  .(p        ] ,    are  a  good  set  to  use  in  a  selective  reduction  treat- 
ment, it  follows  that  we  can  construct  several  reduced  Hamiltonian  ma- 
trices, one  for  each  ip        ,  similar  to  (5.33),  by  simply  taking  the  off- 
diagonal  elements  into  account.   By  so  doing,  we  get  to  order  m 


XX        o    '  o       o      1   '  "^o 


(5.34) 


where 


(i)  _  p  (o,i)   (l,i)    „(?,!). 


(5.35) 


and 


^t(ij)  _  r^(i,j)„(i,j)    rr^^'J^l 


(5.36) 


This  form  corresponds  to  a  partitioned  eigenvalue  problem  in  the  basis 
Ixf    ''     thereby  allowing  complete  flexibility  in  the  determination  of 
the  variation  solution  to  yield  the  best  possible  result.   At  the  same 
time,  the  proper  noninteracting  properties  are  also  guaranteed,  ensuring 
upper  bounds  for  all  roots.   Since  one  may  construct  p  +  1  effective  Ham- 
iltonian matrices  of  the  form  (5.34),  diagonalization  of  each  one  yields 
(p  +  1)  eigenvalues  and  eigenvectors.   From  the  previous  discussion  of 
steepest  descent,  hov7ever,  it  follows  that  by  choosing  i  as  the  particular 
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state  of  interest,  the  i   root  should  be  selectively  reduced  for  m  <  n. 
l«Jhen  m  =  n,  (5.34)  for  any  i  becomes  exactly  equivalent  to  (5.33)  giving 
p  +  1  identical  secular  equations.   Besides  the  rigorous  solutions,  re- 
sults obtained  from  the  approximate  expressions  (5.32)  will  also  be  found 
to  be  useful. 

V.3.  The  Geometric  Formula  for  a  Multi- 
Dimensional  Reference  Space 

In  a  practical  application  of  the  multi-dimensional  reduction  pro- 
cedure, one  is  interested  in  a  CI  problem  in  the  space  J  fE)  partitioned 
into  two  subspaces  js)  and  |tt).   Typically,  js)  is  composed  of  a  selec- 
tion of  n   configurations,  neglecting  the  remaining  n  configurations 
which  are  grouped  into  |  tt ) .   The  dimension  of  |s)  is  usually  chosen  to 
be  much  larger  than  the  number  of  roots  desired.   In  such  a  case,  the 
first  order  solution,  a  type  of  geometric  formula,  has  a  convenient 
simplifying  feature. 

Let  us  consider 

Iff)  =  h  ^)  (5.37) 

for  |ff)an  orthonormal  set  of  configurations.   The  equation  we  would  like 
to  solve  is  (5.4),  or  equivalently, 

e^O;  ^  =  (E  ^E.  (5.38) 

In  practice,  though,  we  usually  must  resort  to  solving  a  truncated  prob- 
lem within  the  chosen  space,  |s).   By  so  doing,  from  the  definitions 
(4.20)  and  (5.37),  we  have 

<s|k|$)  =  <s|h|$>  =  U   ,  (5.39) 


ss 


and 
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HI   ID   =  ID   t  ,  (5.40) 

ss   s      s  o  \  •   / 

for  the  t      of  (5.18)  and  the  n  possible  eigenvectors, 
o  s  * 

^s  =  ^^0^1  •••  ^n  -1^-  (^•^^> 

s 

Since  we  are  only  interested  in  the  first  p  +  1  solutions,  let  us  define 

them  as 

Xi=|Oca^,        i  =  0,l,...p:sng.  (5.42) 

These  wavefunctions  will  constitute  a  multi- dimensional  unperturbed  ref- 
erence space  to  implement  a  perturbation  solution  of  (5.38).  Now  we  may 
partition  |if)  such  that 

Iff)  =  b.^.  t)  =  l^^.ih).  (5.43) 

where  |^)  consists  of  the  remaining  n   -  (p  +  1)  functions,  which  may  be 
obtained  from  the  basis  j$).   The  |y  )  are  assumed  to  be  orthogonal  to 
|>t),  and  an  adequate  realization  would  be  the  other  eigenvectors  of  (5.41), 
which  are  necessarily  noninteracting  with  the  unperturbed  solutions.   The 
follov7ing  is  a  consequence  of  this  latter  property. 

The  first  order  solution  to  the  super-matrix  equation,  (5.22),  is 


XX 


n^:/  =   «^o  +  ^o^s^o  "   'h'^'^'^r  <^5-^^> 


From  the  relation  (5.27),  the  important  components  are  the  [a^  }  de- 
fined  in  (5.26)  where  k  =  0  and  k  =  1.  With  the  definition  of  K  and  by 
using  the  eigenfunction  properties  of  the  |y),  it  follows  that 

.(ij) 


a 
o 


=   <X^|H|lt><tt|H|xj>,  (5.45) 

and 
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o[^'^^    =  <X^lH|tt)<ttlH|u:)<tt|H|x^>,  (5.46) 

for  any  i,j  =  0,1, ...,p.   This  has  the  convenient  consequence  that  in  a 
truncated  CI  calculation  the  first  order  improvement  to  the  p  +  1  lowest 
states  can  be  obtained  by  ignoring  the  n   -  (p  +  1)  combinations,   >(), 
and  solely  introducing  the  excluded  configurations,  |  tt ) ,  into  the  calcu- 
lation.  Since  the  selective  reduction  equations  (5.34)  and  the  approxi- 
mate solutions  (5.32)  are  sub-classes  of  (5.22),  this  feature  is  also 
present  for  these  methods.   Although  we  employed  the  n   -  (p  +  1)  remain- 
ing eigenvectors  of  (5.41)  in  this  derivation,  in  the  interest  of  going 
to  higher  orders,  it  may  be  mentioned  that  this  is  not  required  since 
only  the  orthogonality  of  the  |^  )  to  the  unperturbed  solutions,  |y),  is 
sufficient  to  ensure  the  noninteracting  property. 

\</hen  one  selects  a  set  of  configurations  for  a  CI  calculation,  it  is 
assumed  that  those  remaining  are  relatively  insignificant.   By  determining 
the  first  order  improvement  to  the  initial  choice  |s)  by  means  of  (5.45) 
and  (5.46),  this  supposition  can  be  placed  on  a  more  quantitative  basis. 
Its  assistance  in  choosing  important  configurations  is  also  apparent. 
Furthermore,  since  this  solution  typically  accounts  for  more  than  fifty 
per  cent  of  the  obtainable  energy  improvement,  it  seems  that  it  would  be 
a  logical  extension  of  any  standard  CI  computation. 


CHAPTER  VI.  NUMERICAL  APPLICATIONS  OF  THE 
REDUCED  PARTITIONING  PROCEDURE 


In  order  to  undertake  a  numerical  assessment  of  the  reduced  partition- 
ing theory,  we  made  some  numerical  applications  within  the  configuration 
interaction  framework.   Since  one  of  the  objectives  of  the  theory  is  to 
minimize  the  difficulties  encountered  in  configuration  selection,  we  were 
interested  in  limiting  these  initial  computations  to  molecules  where  the 
full  CI  wavefunction  for  a  basis  of  several  functions  would  be  easily 
tractable,  since  the  full  CI,  being  the  best  answer  in  the  chosen  basis, 
will  provide  the  ultimate  result  to  use  in  our  comparisons.   Practically, 
this  limited  us  to  the  two-electron  systems  H^  and  HeH  ,  as  the  number  of 
initial  basis  functions  allowed  by  our  computer  program  (18)  would  lead 
to  a  few  thousand  configurations  of  the  proper  symmetry  for  even  LiH. 

The  calculations  that  we  have  performed  fall  essentially  into  two 
categories:   ground  state  studies  relative  to  a  single  unperturbed  solu- 
tion, and  a  simultaneous  treatment  of  several  states  by  employing  a  multi- 
dimensional reference  manifold.   The  first  of  these  will  be  illustrated 
with  equilibrium  studies  of  H„  relative  to  a  simple  Coulson  molecular 
orbital  wavefunction  as  the  unperturbed  state  as  well  as  with  a  Hartree- 
Fock-Roothaan  reference  function.   In  addition,  a  complete  potential  curve 

1-4-  4- 

for  the  lowest   Z   state  of  HeH   subject  to  a  liFR  unperturbed  solution 

will  be  obtained.   As  an  example  of  the  multi-dimensional  case,  the  selec- 

1  +  + 

tive  reduction  of  the   E  excited  states  for  HeH  v^ill  also  be  considered, 

including  some  numerical  justification  of  the  approximate  treatment  pre- 
viously described.   Although  this  discussion  is  self-contained  in  the 
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sense  that  the  objective  is  a  theoretical  solution,  due  to  the  consider- 
able experimental  interest  in  the  reactions  of  metastable  He  [55,56], 
these  calculations  also  may  be  of  some  interest  in  connection  with  this 
problem. 

VI. 1.  Summary  of  the  Computational  Procedure 
The  computer  program  implemented  for  this  problem  couples  the  in- 
gredients of  a  configuration  interaction  calculation  with  the  perturba- 
tion approach  described  in  the  previous  chapters.   The  initial  basis 
functions  used  are  the  complex  form  of  the  Slater  type  orbitals  (STO's) 
located  on  each  center  in  the  molecule.   The  complex  form  ensures  that 
each  orbital  is  an  angular  momentum  eigenstate  and  thus  assists  one  in 
constructing  n-particle  functions  of  the  proper  s>'mmetry.   From  the 
initial  basis  of  STO's,  a  secondary  basis  consisting  of  an  orthogonal  set 
of  orbitals  is  generated.   In  the  calculations  presented  here,  these 
are  either  the  Schmidt  orthogonalized  combinations  of  the  STO's  or  the 
Hartree-Fock-Roothaan  eigenvectors  one  obtains  from  solving  a  preliminary 
single  determinant  problem.   By  associating  a  or  |3  spin  with  each  member 
of  this  secondary  basis,  a  set  of  one-particle  pure  spin  orbitals  is  ob- 
tained.  Due  to  the  symmetry  of  the  STO's,  the  resultant  spin  orbitals 
are  specified  by  the  usual  a,  rr,  6,  etc.  classification.   From  the  spin 
orbitals,  n-particle  anti-symmetrized  products  are  constructed,  where  the 
symmetry  of  the  state  S,  11,  A,  etc.  is  easily  achieved  from  the  spin  orbi- 
tal specifications.   For  homonuclear  diatomics  the  gerade  and  ungerade 
symmetry  for  the  orbitals  is  also  taken  into  account.   By  combining  these 
anti-symmetrized  products  into  the  proper  combinations,  spin  symmetry  is 
incorporated  into  the  elements  of  the  basis  set,  as  is  the  plus  and  minus 
sjTiimetry  arising  from  reflection  in  the  plane  of  the  molecule  for  sigma 
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states.   A  configuration  for  a  diatomic  molecule  is  thus  defined  as  a 
pure  s>'Tiimetry  function  in  all  respects,  as  this  gives  the  minimum  possi- 
ble number  that  must  be  considered  for  the  problem  of  interest. 

The  molecular  integrals  are  computed  from  a  program  written  by  H.  H. 
Michels  and  altered  by  J.  Kouba  to  give  integrals  accurate  to  at  least 
eight  figures.   Additional  routines  from  J.  Kouba' s  CI  program  [57]  have 
been  modified  by  the  author  for  inclusion  into  the  completed  program. 
I'Jhen  Hartree-Fock  wavefunctions  are  required,  they  are  obtained  from  a 
program  written  by  H.  H.  Michels,  but  with  the  more  accurate  integral 
routines  substituted. 

It  may  be  worthwhile  to  briefly  summarize  the  essential  steps  in  the 
remainder  of  the  program,  which  is  used  to  effect  the  perturbation  calcu- 
lations.  This  section  is  arranged  to  calculate  CI  wavefunctions  relative 
to  a  given  set  of  configurations,  and  then  make  two  sub-partitionings  in 
each  single  run.   With  each  partitioned  result  constituting  one  or  several 
-unperturbed  states,  these  solutions  and  the  remaining  elements  of  the 
orthogonal  complement  are  used  to  construct  the  perturbation  quantities 
by  the  theory  previously  described.   In  this  respect,  routines  have  been 
written  which  allow  us  to  work  either  within  the  RS  or  BW  framework. 
Since  a  variational  solution  relative  to  the  perturbation  basis  is  desired, 
this  can  be  accomplished  by  generating  the  respective  Hamiltonian  and 
overlap  matrices  from  which  the  eigenvalue  equation  may  be  solved.   In  the 
BW  case,  this  is  easily  accomplished  solely  from  the  perturbation  energies, 
since  from  the  basic  definitions,  (2,40),  it  follows  that 


and 


<cp    Icp   )  -   e        .  (6,1b) 

p    q        i>i-q  V  •    -» 
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In  the  RS  case,  the  overlap  matrix  elements  must  be  calculated  from  the 
perturbation  corrections  (2.52),  but  from  these  overlap  elements  and 
the  perturbation  energies  (2,53),  the  Hamiltonian  matrix  elements  are 
efficiently  constructed  via  a  relationship  derived  by  Lowdin  [l5(PTTX)] 
giving 

<$  V  $  )  -  G    ,  +  Z       S   €    ,  ,  ,<$,  U,).     (6.2) 

Within  BW  theory,  an  iterative  procedure  for  obtaining  the  [M,M-1] 
and  [m,M]  Fade'  approximants  by  means  of  the  solution  of  the  implicit 
equation,  £■  =  f(S),  has  also  been  implemented.   This  routine  incorpoTates 
the  technique  of  Newton-Raphson  iterations  [15(PTI)],  a  second  order  pro- 
cess, to  increase  the  rate  of  convergence.   For  example,  for  the  [m,M-1] 
approximant  with  a  single  partitioning 

e  =  f(£)  =  E°  +   c?''*(eE  -  E,  )""^(2r  ,  (6.3) 

ooolo 

from  which  one  readily  obtains  the  quantity 

f'(£)  =  -(^''"(eE   -  E^)"-^E  (£E   -  E.y'^d     ,     (6.4) 

GO       1        O      O       i       O 

which  is  always  negative.   From  a  first  order  approximation 

fi^'^  =  f(£^°')  .  (6.5) 

the  improved  value, 

&*  =  £(°)  +  [(£^1)  .  e(°>)/(l  -  f(£^°>)],  (6.6) 

follows  from  an  application  of  the  Newton-Raphson  formula  to 

F(£)  =  £  -  f(£)  =  0.  (6.7) 

Furthermore,  the  fact  that  f '  (£)  <  0  guarantees  that  for  S,^'^^    >  £, 
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(I)  * 

£^  '  <  5;  but  S  ,  v/hich  can  be  ralatcd  to  the  variational  principle 

[15(PTI)],  must  always  be  greater  than  £,  ensuring  convergence  from 

above.   This  process  is  continued  until  two  successive  £' s  differ  by  no 
,-10 


mo 


re  than  10    .   Since  this  method  requires  repeated  inversions  of 


(£  E  -E^)  for  the  various  values  of  £,  although  it  may  be  more  desirable 
o   1 

in  some  respects,  it  is  a  more  time  consuming  process  than  the  equivalent 
secular  equation  solution. 

Despite  the  particular  approach  used,  one  encounters  computational 
problems  V7ith  near  linear  dependency  in  the  perturbation  basis  set.   In 
the  implicit  function  process,  tliis  complication  is  manifested  in  the 
inversion  of  the  ill-conditioned  matrix  (SE  -E.),  and,  in  the  eigenvalue 
equation,  by  the  removal  of  the  overlap  matrix.   In  the  program,  this 
possibility  is  investigated  at  each  order  of  calculation  by  multiplying 
the  matrix  with  its  inverse,  or  checking  the  traiisformed  overlap  matrix. 
As  the  perturbation  functions  rapidly  exhaust  the  space  of  interest,  a 
resulting  near  linear  dependency  is  indicative  of  having  essentially 
achieved  the  desired  solution,  but  in  the  process   it  is  often  of  interest 
to  gain  another  one  or  two  figures  of  accuracy.   Consequently,  to  obtain 
slightly  higher  orders  of  solution  than  v/ould  otherwise  be  possible,  a 
multi-precision  extension  written  by  C .  E.  Reid  [58]  for  the  IBM  360/65 
has  been  iiicorporated  into  the  program.   With  this  package,    multi- 
precision  routines  for  the  Fade  approxiraant  analysis,  including  a  matrix 
inversion  subroutine,  have  been  written.   The  use  of  these  routines  de- 
pends on  the  assumption  that  the  pertinent  quantities  need  to  be  calcu- 
lated only  in  double  precision  (16  figures).   Then  with  these  values,  the 
remaining  analysis  involving  the.  more  sensitive  steps  is  accomplished 
with  four  words  (37  figures)  of  precision.   This  process  can  be  trusted 
up  to  the  point  where  the  linear  dependency  once  again  manifests  itself. 
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or  when  the  ultimate  effect  depends  on  differences  beyond  the  range  of 

the  low  precision  of  the  initially  calculated  quantities.   Either  of 
these  features  becomes  apparent  in  the  course  of  the  calculation. 

In  order  to  treat  multiple  states  with  the  selective  reduction  pro- 
cedure, the  program  uses  essentially  the  same  computational  framework 
as  in  the  single  partitioning  case.   To  illustrate  this,  let  us  consider 
the  secular  equation  approach  which  requires  constructing  the  "reduced" 
Hamiltonian  and  overlap  matrices. 


HI 


(i) 
R 


and 


(i) 


The  partitionings  are  defined  as 


HI 


HI 


H 


XX 

(i) 

X9 

(i) 

cpcp 


=     t 


,(i) 


(i,0)„(i,l) 


m  -' ' 


(6.8) 


(6.9) 


(6.10a) 
(6.10b) 


and 


^(i)        ,    (i)|„(i) 


-{'r^^'^) 


(i) 


(6.11a) 
(6.11b) 
(6.11c) 


for  the  reference  functions,  [y),and  the  sigma  terms  defined  in  (4,49), 
(4.50),  and  (5.26), 
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Since  we  have  t       from  the  unperturbed  solution,  the  quantities 


o 


needed  are  HI    ,  H    ,  and  S.         .      To  obtain  these  matrices,  it  is  con- 
venient  to  begin  by  calculating  the  approximate  solution,  (5.32),  for 
each  of  the  p+  1  states.   As  this  is  a  single  partitioning  approach,  this 
is  achieved  with  the  basic  program  structure.   In  particular,  the 
quantities 

b.  =  <x.|K|lh)  (6.12) 

and  the  perturbation  functions 

p^'^  =  |ll.>l3^'^  (6.13) 

are  generated  for  each  i.   From  these,  one  also  produces  the  matrices 

H     and  S.         ,  v/hich  defines  the  secular  equation  for  the  approximate 

cpcp  cprp  ' 

solution,  (5.32).   By  solving  this,  the  approximate  energy,  S.  ,  and 
wavefunction,  ijf.,  are  given  as  the  lowest  energy  solution.   [n  the  pro- 
cess of  obtaining  each  of  these  approximate  results,  by  simply  storing 

the  lj3    ,  lb  .  ,  HI    ,  and  A    externally,  all  the  basic  information 
1-    W  '       cpcp 

necessary  for  the  rigorous  multiple  states  solutions  is  available. 

The  remaining  elements  required  are  the  terms  {o     }  which  arise 
from  the  interaction  of  the  unperturbed  states,  and  are  thus  unique  to 
the  multi-dimensional  case.   These  present  no  problem,  though,  since  from 
the  definition 

^k-1^^  =  <X.|K|cpj^'^,      k  ^  1,  (6.14) 

one  immediately  obtains 

^k-l^^  =  (b  .  Ii3^^^  .  (6.15) 

Following  this  strategy,  the  program  efficiently  constructs  the  p  +  1 
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possible  reduced  Hamiltonian  and  overlap  matrices,  from  vjhich  the  secular 
solution  gives  p  +  1  roots,  of  which  the  i   one  is  selectively  redxiced. 
As  a  bonus  the  approximate  solutions  are  given,  and  as  our  computations 
will  indicate,  these  are  likely  to  be  sufficiently  close  to  constitute 
a  satisfactory  treatment  for  many  applications. 

To  place  the  computational  benefits  of  the  reduced  partitioning 
theory  into  perspective,  it  should  be  noted  that  in  any  truncated  CI  cal- 
culation with  a  basis  set,  the  calculation  of  the  molecular  integrals  Is 
the  predominant  time  consuming  operation,  with  a  considerable  amount  of 
time  also  required  for  their  transformation  to  a  secondary  basis.   Sina?.  .  ■ 
both  of  these  steps  are  necessary  despite  the  number  of  configurations  •- 
used  in  the  calculation,  it  seems  logical  to  get  maximum  benefit  from  the 
chosen  basis  by  using  many  more  configurations  than  is  customary,  or  evea 
all  that  are  possible.   To  do  this  in  the  standard  approach  would  involve 
solving  an  N  X  N  secular  equation  for  several  states,  v^7hich  is,  at  pres- 
ent, an  impossibility  for  very  large  N.   On  the  other  hand,  the  reduction 
process  does  not  require  such  a  diagonalization,  but,  as  will  be  shown, 
gives  a  result  in  essential  agreement  with  the  full  solution  by  only 
solving  low  order  secular  equations.   It  is  still  generally  necessary  Lu 
construct  all  the  Hamiltonian  matrix  elements,  but  even  though  this  is 
admittedly  a  problem  for  very  large  N,  since  these  elements  are  obtainable 
solely  from  the  transformed  molecular  integrals  and  the  configurations, 
the  extension  is  reasonable. 

* 

VI. 2.  Results  and  Discussion:   The  Single  Partitioning  Case 

As  a  first  example  of  the  reduced  partitioning  procedure  (RPP) ,  we 
considered  the  hydrogen  molecule  at  the  equilibrium  internuclear  separation 


"All  the  calculations  reported  herein  were  done  on  the  IBM  360/65 
computer  of  the  University  of  Florida  Computer  Center. 
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R  =  1,40  a.u.   The  n-particle  basis,  |cf)  Ih),  V7as  generated  from  an  ini- 
tial set  of  nine  STO's  on  each  atom,  ls(1.20),  l5'(1.00),  23(1.16), 
2p  (1.71),  2p   (1.71),  3d  (2.20),  and  3d   (2.00),"  which,  using  Schmidt 
orthogonalized  combinations,  led  to  34  possible  configurations  of   E 
symmetry.   As  this  is  the  full  CI  basis,  it  gives  the  best  answer  obtain- 
able for  the  chosen  space.   The  full  CI  result  is  E  =  -1.1712787  a.u. 
The  principal  contribution  to  the  CI  wavefunction  is  the  Coulson 
function  configuration 

^n  V'Al^^^^l)  Is  0(2)  )  (6.16) 

which  constitutes  ■«  987o  of  the  full  solution  and  gives  E  =  -1.1280771  a.u. 
by  itself.   Choosing  this  configuration  as  the  reference  function,  the 
[cp  }  and  the  [o.]   were  calculated  from  (4.44)  and  (4.45),  and  from  these 

I^  J. 

quantities  the  reduction  energies  were  obtained  from  (4.51)  for  each  m 
from  1  to  13.   These  numbers  are  listed  in  Table  1. 

For  comparison,  the  difference  between  the  unperturbed  energy  and 
the  full  CI  result  will  be  defined  as 

AE(i)  =  e°  -  E.  (6.17) 

where  i  indicates  the  state  of  interest.   If  £.  is  a  HF  energy,  then  AE 

may  be  referred  to  as  E     ,  the  correlation  energy  obtainable  within  the 
■'  corr 

space  |ff),  by  analogy  with  the  usual  definition  of  correlation  energy. 

The  most  dramatic  feature  illustrated  in  Table  1  is  the  anticipated 
extreme  divergence  of  the  sigma  series  which,  lacking  a  non- linear  summa- 
tion technique,  would  be  useless.   However,  by  applying  the  [M,M-1]  Fade 


"The  orbital  exponents  for  n  =  2  were  taken  from  the  optimized  basis 
set  for  R  =  1.40  a.u.  obtained  by  McLean,  £t  al.  [59], 

•'■"In  fact,  Coulson  [60]found  the  best  scale  factor  equal  to  1.197 
for  R  =  0.732  A. 
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TABLE  1 

Reduction  Energies  for  H„  with  the  Coulson  Reference  Function" 
-^  (R  =  1.40) 


,(m) 


M     £^'"^  =  e^+  l/„(m)[M,M-l]    %  e; 

o     o  coir 


Brillouin-Wigner 
Sigma  Elements 


0 

-1.1230771 

0.00 

^0  = 

0.124303 

1 

-1.1582136 

69.76 

^1  = 

0.279955 

2 

-1.1662436 

88.35 

^2  = 

1.594896 

3 

-1,1696479 

96.23 

-3- 

1.25  X  10^ 

4 

-1.1707700 

98.82 

^4  = 

1.18  X  10^ 

5 

-1.1711049 

99.60 

^5  = 

1.22  X  10^ 

6 

-1.1712244 

99.87 

• 
• 

7 

-1.1712552 

99.95 

^25 

23 
-  2.23  X  10 

8 

-1.1712685 

99.98 

9 

-1.1712741 

99.99 

10 

-1.1712764 

99.99 

11 

-1.1712778 

100.00 

12 

-1.1712782 

100.00 

13 

-1,1712784 

100.00 

33(Full  CI) 

-1.1712787 

Energies  and  bond  distances  are  always  expressed  in  atomic  units, 
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approximants  to  the  series,  it  is  seen  that  the  reduction  energies  con- 
verge very  rapidly  toward  the  full  CI  solution  [61].   From  the  connection 
between  the  [M,M-1]  approximant  and  the  variational  treatment,  it  is  more ■ 
easily  seen  why  this  result  is  found  [61].   This  is  the  same  situation 
previously  observed  in  an  application  of  BW  perturbation  theory  to  the 
hydride  ion  [52],  where  Pad^  approximants  were  also  used  to  sum  the  sericj 
to  give  meaningful  energies. 

As  has  been  mentioned,  the  fact  that  the  perturbation  functions 
[cp  ,cp,  , .  .  .  ,cp  )  rapidly  exhaust  the  space  of  interest  causes  problems 
with  near  linear  dependency.   In  the  examples  studied,  this  usually  be- 
comes apparent  in  the  ninth  or  tenth  order  solution.   Consequently,  mul- 
tiple precision  was  used  to  obtain  the  higher  order  solutions  in  Table  1, 
permitting  another  one  or  two  digits  of  agreement. 

Due  to  the  simplification  inherent  in  the  first  order  solution,  it 
is  especially  encouraging  that  this  result  accounts  for  almost  seventy 
per  cent  of  the  possible  energy  improvement.   For  many  problems  it  may  be 
difficult  to  calculate  the  higher  order  terms,  thus  a  dramatic  energy  de- 
crease from  the  "geometric"  result  is  a  definite  asset. 

In  order  to  consider  a  IIFR  reference  function,  the  previous  STO  basis 
set  was  used  in  a  IIFR  calculation  to  give  a  set  of  one-electron  orbitals 
as  eigenfunctions  to  the  Fock  Hamiltonian.   From  these,  another  34-config- 
uration  CI  wavefunction  was  constructed  but  with  the  principal  configura- 
tion now  actually  being  the  HF  result.   The  energy  of  the  HF  solution  is 
-1.1335224  a.u.,  and,  of  course,  the  full  CI  energy  is  the  same  as  before. 

For  a  molecule  such  as  H„ ,  where  the  Is   STO  combination  is  a  good  approx- 

2'  g 

imation  to  the  1^  molecular  orbital,  there  cannot  really  be  too  much  dif- 
S 

ference  between  the  Coulson  function  and  the  HF  solution.   This  is  reflect- 
ed in  a  comparison  of  Table  1  and  Table  2.   The  elements  of  the  sigma 
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TABLE  2 

Reduction  Energies  for  H»  with  a  Hartree-Fock  Reference  Function 

(R  =  1.40) 


M     e*"™^  =  e  +  l/.(m)[M,M-l]     %  E^ 

o     £.^  ^'-  '    -'       corr 


Brillouin-Wigner 
Sigma  Elements 


0 

-1.1335224 

0.00 

^0  = 

0,106462 

1 

-1.1590861 

67.71 

^1  = 

0.243929 

2 

-1.1668000 

88.14 

^2  = 

1.509162 

3 

-1.1696330 

95.64 

^^3  = 

1,26  X  10^ 

4 

-1.1707910 

98,71 

^4  = 

1.22  X  10^ 

5 

-1.1711357 

99.62 

°5  = 

1.28  X  10^ 

6 

-1.1712490 

99.92 

• 
t 
• 

7 

-1.1712697 

99.98 

^23 

22 

=  1.27  X  10 

8 

-1.1712763 

99.99 

9 

-1.1712779 

100.00 

10 

-1.1712785 

100.00 

11 

-1.1712786 

100.00 

12 

-1.1712786 

100.00 

33(Full  CI) 

-1.1712787 

74 

expansion  are  very  similar  as  are  the  reduction  energies  for  each  order, 
although  the  results  with  the  HF  reference  function  do  seem  to  converge 
a  little  better  in  higher  order.   The  energy  decrease  given  by  the  geo- 
metric solution  is  somewhat  better  for  the  Coulson  function,  which  is 
probably  a  result  of  more  improvement  being  possible.   It  should  be  ob- 
served that  better  than  ninety-five  per  cent  of  the  possible  correlation 
energy  is  given  by  only  a  third  order  or  four-function  expansion, 

For  the  HeH  molecular  ion  a  very  good  HF  calculation,  including 
orbital  exponent  optimization,  has  been  done  by  Peyerimhoff  [63].   For 

R  =  1.455  a.u.,  the  optimized  basis  she  obtained  is  ls„  (1.37643), 
'       t-  He 

Is'   (3.87107),  2s,^  (1.54335),  2p  ^  (2.64576),  2p'  ,^  (3.24082),  3d  .^ 
He  He  o  He  ^o  He  o  He. 

(2.54147),  4f  „  (3.73526),  Is,  ri .00949) ,  2s„(l . 18036) ,  2s '  (2.56229) , 
o  He  n  H  H 

2p   ,,(1.79089),  3d  ,,(2.41228),  from  which,  she  reports  an  energy  of  E,^_  = 
o  H  o  H  Hr 

-2.933126  a.u.   By  repeating  this  HF  calculation,  v/e  obtained  the  slightly 

higher  result,  E   =  -2.933072  a.u.,  which  we  expect  is  due  to  our  more 
HF 

accurate  integral  computation.   By  augmenting  this  sigma  basis  with  2p 

±1  He 

(2.868),  2p  ,  ,,(0.827),  and  3d  ^  „  (2.686),"  a  full  CI  wavefunction  con- 
'^  ±1  H  ±2  He 

1  -4- 
sisting  of  82  configurations  of  E  symmetry  was  obtained.   With  the  HF 

solutions  as  the  reference  function,  the  series  elements  and  the  reduction 
energies  are  as  given  in  Table  3.   It  appears  the  sigma  series  is  even 
somewhat  more  divergent  than  in  the  H-  case  ^^7hich  might  be  expected  from 
the  convergence  criterion  given  in  (IV. 3),  but  may  also  be  dependent  on 
the  adaptibility  of  the  basis  set,  or  the  number  of  configurations.   In 
keeping  with  the  greater  extent  of  divergence,  the  convergence  of  the  re- 
duction energies  is  slightly  worse  than  before,  even  though  the  third 
order  solution  still  accounts  for  more  than  ninety  per  cent  of  the  available 
correlation  energy.   The  geometric  result  gives  about  fifty-nine  per  cent. 


"■This  STO  basis  will  be  referred  to  as  basis  I, 
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TABLE  3 

Reduction  Energies  for  HeH  with  a  Hartrec-Fock  Reference  Function 

(R  =  1.455) 


M 


S^""^  =  £  +  l/„(m)[M,M-l]    %  E^ 

o     S^  ^'-  '    -'       corr 


Brillouin-Wigner 
Sigma  Elements 


0 

-2.933072 

0.00 

^0=-- 

0.260175 

1 

-2.956818 

59.21 

'^l  = 

1.723745 

2 

-2.966665 

83.77 

^2  = 

3.01  X  10^ 

3 

-2.969824 

91.64 

^3  = 

6.58  X  10^- 

4 

-2.971765 

96.48 

^4  = 

1.64  X  10^ 

5 

-2.972594 

98.55 

^5  = 

4.39  X  10^ 

6 

-2.972941 

99.41 

• 
• 
* 

7 

-2.973046 

99.68 

^23 

=  5.60  X  10^° 

8 

-2.973105 

99.82 

9 

-2.973138 

99.90 

10 

-2.973153 

99.94 

11 

-2.973162 

99.97 

12 

-2.973168 

99.98 

81  (Full  CI) 

-2.973176 
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To  obtain  a  more  balanced  basis  set  for  potential  curve  calcula- 


-f. 

tions  for  HeH  ,  basis  II  was  obtained  by  excluding  the  tt  and  6  functions, 

as  well  as  the  4f  (lie)  and  2s '  (H)  functions,  ivhich  are  least  important 
to  the  HF  solution.   These  were  replaced  by  2p      (3.00),.  2p     (2.00), 

3d  2  u  (3.00),  and  3d  „   (2.00).   This  basis  results  in  61  configurations 

1  -f- 

of     E      symmetry.      The  HF   energy  at   R  =   1.460   a.u.    is    -2.93302   a.u.      In 

1  j- 

Figure  1,  the  HF  potential  curve  and  the  full  CI  curve  for  the  X  E 

ground  state  can  be  compared.   Since  He  is  a  closed  shell  atom,  this  is 
one  of  the  few  special  cases  when  the  HF  potential  separates  properly, 
giving  a  HF  atom  and  a  proton.   The  main  feature  to  be  observed  is  the 
relative  constancy  of  the  correlation  correction.   In  Table  4  the  reduc- 
tion energies  through  the  eighth  order  are  given  as  a  function  of  R. 
From  the  percentages  of  the  correlation  energy  listed,  the  convergence 
with  only  nine  functions  is  seen  to  be  quite  satisfactory,  being  poorest 
at  R  =  2.50  a.u.  and  excellent  at  large  internuclear  separations.   The 
first  order  solution  accounts  for  better  than  seventy  per  cent  for  larger 
R  values.   Figure  2  illustrates  this  as  well  as  the  additional  improvement 
given  by  the  second  order  solution.   It  should  be  noted  also  that  the  cor- 
relation energy  is  not  actually  constant  but  shows  a  rise  at  R  =  1.52  a.u. 
and  decreases  slightly  as  R  increases. 

The  question  of  the  shapes  of  the  potential  curves  for  the  reduction 
energies  may  also  be  of  importance  if  one  is  interested  in  the  vibrational 
spectra.   Thus,  to  more  adequately  assess  the  differences  betvjeen  the 
curves  in  Figure  2,  a  Dunham  analysis  was  applied  to  each  of  these  approx- 
imations using  a  program  written  by  Beebe  [64].   These  results  are  pre- 
sented in  Table  5,  along  with  some  comparison  values  obtained  in  other 
calculations  by  Peyerimhoff  [63],  Anex  [65],  Michels  [66],  and  Wolniewicz 
[67],   Since  it  is  not  necessarily  obvious  just  how  to  treat  the  separated 


77 


,(m) 


TABLE  4 

1  -)-  + 

Reduction  Energies  for  the  X  E  State  of  HeH 

as  a  Function  of  Internucleai"  Separation^ 
(Single  partitioning  with  a  HF  reference  function) 


1.00 


1.40 


1.46 


1.52 


^0 

2.86018 

2.93249 

2.93302 

2.93240 

e^^) 

2.88522 

2.95771 

2.95820 

2.95733 

e(^) 

2.89392 

2.96703 

2.96757 

2.96693 

e(3> 

2.89705 

2.97071 

2.97130 

2.97070 

£(^> 

2.89789 

2.97179 

2.97243 

2.97188 

£(5) 

2.89807 

2.97205 

2.97270 

2.97217 

e(^> 

2.89814 

2.97218 

2.97283 

2.97231 

e(^) 

2.89817 

2.97223 

2.97289 

2.97237 

&(^) 

2.89818 

2.97226 

2.97292 

2.97240 

E  (Full  CI)    2.89818 
o 


2.97228      2.97295 


All  energies  are  negative. 


2.97243 


AE 

0.03800 

0.03979 

0.03993 

0.04003 

7„AE(e*^^^ 

65.89 

63.38 

63.06 

62.78 

7oAE(e^^^) 

100.00 

99.95 

99.92 

99.93 
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.(m) 


TABLE  4  (continued) 


2,00 


2.50 


3.50 


5.00 


«0 

2.91057 

2.88776 

2.86746 

2.86218 

e<'' 

2.93533 

2.91243 

2.89242 

2.88739 

e<2> 

2.94456 

2.92112 

2.90039 

2.89507 

e"> 

2,94835 

2.92454 

2.90304 

2.89752 

e<*> 

2.94964 

2.92557 

2.90371 

2.89808 

e<5> 

2.94994 

2.92578 

2.90387 

2.89819 

£<^> 

2.95011 

2.92592 

2.90393 

2.89822 

e<'>' 

2.95021 

2.92600 

2,90395 

2.89822 

e<'' 

2.95026 

2.92604 

2.90396 

2.89823 

E  (Full  CI)   2.95033 
o 


2.92612      2.90398 


2.89823 


AE 

0 

.03976 

0.03836 

0,03652 

0,03605 

7oAE(£^^^) 

62.27 

64.31 

68.35 

69.93 

7=AE(e^^^) 

99.82 

99.79 

99.95 

100.00 
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TABLE  4  (continued) 


,(ra) 


6.00 


8.00 


10.00 


«0 

2.86169 

2.86147 

2.86142 

£<^' 

2,88694 

2.88673 

2.88668 

£(2> 

2.89458 

2.89435 

2.89430 

e(^> 

2.89702 

2.89680 

2.89678 

£("> 

2.89756 

2.89733 

2,89727 

£<5) 

2,89767 

2.89743 

2.89738 

fiC^' 

2.89769 

2.89746 

2,89740 

£<') 

2.89769 

2.89746 

2.89741 

8<«) 

2,89770 

2.89746 

2.89741 

E  (Full  CI)   2.89770 
o 


2,89746      2.89741 


AE  0,03601 

%AE(e*-^^)      70,12 
%AE(e^^^)     100.00 


0.03599 

0.03599 

70.19 

70.19 

100.00 

100.00 
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atom  limit  for  the  £    and  £    curves,  the  values  calculated  at  10  a.u. 
V7ere  assumed  to  be  a  sufficiently  close  estimate  to  use  in  the  deter- 
mination of  the  dissociation  energies.   Anex  and  Wolniewicz  used  the 
exact  value  for  the  He  atom,  E  =  -2.90372  a.u.,  while  in  Michels'  study, 
the  calculated  value  of  E  =  -2.87574  a.u.  was  the  reference  point.   The 
values  for  10  a.u.  should  be  very  close  to  the  calculated  atom  limit, 
but  even  though  our  dissociation  energies  seem  to  be  somewhat  greater 
than  the  comparison  values,  if  any  error  is  present,  the  magnitude  should 
be  slightly  underestimated.   This  is  not  unexpected,  however,  since  bind- 
ing energies  defined  in  this  manner  are  not  governed  by  the  variational 
theorem.   The  lowest  total  energy  found  is  that  of  !\'olniewicz  where  E  = 
-2.97867  at  the  minimum  given  in  Table  5,   Anex's  result  is  E  =  -2.97424, 
while  Michels  reports  a  value  of  E  =  -2.94373  a.u. 

From  the  previous  spectroscopic  parameters  in  Table  5,  it  is  indi- 
cated that  despite  the  dramatic  energy  decrease,  the  actual  shapes  of  the 
potential  curves  for  the  first  and  second  order  reduction  energies  appear 
to  be  somewhat  further  from  the  "best"  results  than  is  the  HF  curve.   This 
effect  of  some  energy  improvement  destroying  agreement  with  other  proper- 
ties is  not  uncommon  in  varieties  of  perturbation  and  CI  studies.   In  addi- 
tion to  the  shape  of  the  potential  curves,  another  manifestation  of  this 
lies  in  the  accuracy  of  the  density  distribution  which  is  usually  apparent 
from  predictions  of  the  dipole  moment  or  other  one-electron  properties. 

Although  the  IIF  wavefunction  is  knovm  to  have  certain  stability  fac- 
tors [8],  for  more  general  cases  these  questions  have  not  been  adequately 
pursued.   However,  one  can  typically  expect  better  agreement  for  proper- 
ties more  sensitive  than  die  energy  when  additional  flexibility  is  built  into 
the  trial  wavefunction.   Since  these  low  order  reduction  solutions  have 
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been  obtained  with  only  one  and  two  variable  parameters,  the  variational 
principle  apparently  overly  prejudices  the  calculation  toward  an  energy 
reductioii  with  a  consequent  sacrifice  in  the  accuracy  of  the  potential 
curves.   To  avoid  these  problems,  somewhat  higher  orders  must  be  included. 
As  these  are  incorporated,  the  near  equality  of  the  eighth  order  sol.itioa 
and  the  full  CI  necessitates  that  the  spectroscopic  parameters  be  essen- 
tially identical. 

In  all  of  the  examples  presented  in  this  section,  a  dramatic  energy 
decrease  was  given  by  the  low  order  solutions   and  essential  agreement 
in  the  higher  orders.   In  each  case,  the  reduction  objective  was  met  with 
8  to  12  functions,  even  though  the  number  of  configurations  varied  from 
34  to  82,   From  the  nature  of  the  procedure,  it  is  expected  that  a  similar 
number  of  functions  should  be  adequate  even  with  many  more  configurations. 
Since  we  can  have  only  single  and  double  excitations  in  the  systems 
studied,  and  the  first  order  solution  relative  to  a  IIF  reference  function 
depends  only  on  the  latter,  it  may  be  argued  that  the  energy  reduction 
obtained  from  the  geometric  formula  may  not  be  as  pronounced  for  larger 
molecules.   However,  since  it  is  generally  conceded  that  from  an  energy 
criterion  the  double  excitations  are  the  most  significant,  we  still  an- 
ticipate a  substantial  improvement  over  the  HFR  result  even  for  other 
molecules.   To  obtain  reliable  spectroscopic  properties,  though,  it  is 
important  to  obtain  slightly  higher  order  solutions.   This  feature  should 
also  apply  to  observables  that  are  determined  by  averaging  over  the 
electron  density. 
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VI  .3.  Results  and  Discussion:   The  Selective  Reducl:lon 
Procedure  for  Multiple  States 

In  order  to  examine  the  selective  reduction  treatment  for  excited 
states,  the  previous  61  configuration,  basis  II  calculation  was  parti- 
tioned to  give  a  ten- function  reference  space  composed  of  the  HF  solution 
and  the  nine  possible  singly  excited  configurations.   By  solving  a  CI 
problem  with  these  ten  functions,  a  set  of  unperturbed  solutions,  jy)  , 
were  obtained.   From  these  reference  functions,  the  selective  process  of 
(V.2)  was  used  to  calculate  reduction  energies  for  the  ground  state  and 
the  four  lowest  excited  states  of  E  symmetry. 

The  selectivity  of  the  reduction  technique  is  demonstrated  in  Table 
6.   Consistent  with  the  steepest  descent  argument,  the  percentages  of 
AE(i)  indicate  that  for  each  choice  of  unperturbed  solution  only  the  cor- 
responding root  is  substantially  affected.   The  most  efficient  reduction 
is  found  in  the  ground  state,  although  the  reduction  process  for  the  low- 
est three  excited  states  is  also  effective.   In  each  of  these  cases,  the 
first  order  solution  is  still  seen  to  be  significant. 

The  ground  state  perturbation  solution  obtained  with  the  multi- 
dimensional reference  space  is  interesting,  since, by  Brillouin's  theorem, 
the  single  excitations  cannot  mix  with  the  HFR  ground  state,  iu-dleave  one 
with  exactly  the  same  unperturbed  energy  as  in  the  single  partitioning 
case.   Since  more  functions  are  explicitly  employed  in  the  multi-dimensional 
approach,  the  increased  flexibility  must  of  necessity  yield  some  improve- 
ment over  the  single  partitioning  case  for  each  order.   Hovrevcr,  the  energy 
decrease  is  found  to  be  negligibly  limited  to  the  fifth  decimal  place,  thus 
the  results  in  Tables  4  and  5  are  the  same  for  all  practical  purposes. 
This  consequence  attests  to  the  fact  that  the  extent  of  interaction  among 
the  elements  of  the  reference  space  and  the  contribution  from  the 


TABLE  6 


,+  a 


Illustration  of  Selective  Reduction  for  Multiple  States  of  HeH 
(Ten-function  partitioning;  HF  and  single  excitations;  R  =  1,460) 
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0 

1 

2 

3 

4 

Unperturbed  -„o 
Energy        1"^ 

2.93302 

1.91869 

1.41766 

1.10141 

0.61147 

Reference 

Function 

0 

>^0 

&(^) 

2.95830 

1.91938 

1.41900 

1.10398 

0.61i9i 

7oAE(i) 

63.31 

00.83 

02.13 

04.79 

00.27 

^(8) 

2.97294 

1.96910 

1.41954 

1.12804 

0.69951 

7oAE  (i) 

99.97 

60.95 

02.99 

49.62 

53.52 

o 

^1 

e(^) 

2.93445 

1.95689 

1.41845 

1.10753 

0.61159 

%AE (1) 

03.58 

46.19 

01.26 

11.40 

00.07 

e(«) 

2.95028 

2,00131 

1.42300 

1.11040 

0.73366 

7oAE(i) 

43.22 

99.89 

08.50 

16.75 

74.28 

o 
^2 

£<^) 

2.93651 

1,91977 

1.44910 

1.10166 

0.61164 

7oAE(i) 

08.74 

01.31 

50.06 

00.47 

00.10 

^(8) 

2,93887 

1.92574 

1.48031 

1.11245 

0.62182 

7=AE(i) 

14.65 

08.52 

99.76 

20.57 

06.29 

o 

^3 

e(^) 

2.93630 

1.92222 

1.41778 

1.12206 

0.61147 

7oAE(i) 

08,21 

04.24 

00.19 

38.48 

00,00 

e^«> 

2.95176 

1,93000 

1.43357 

1,15461 

0.62356 

7oAE(i) 

46.93 

13,67 

25.33 

99.12 

07.35 

All  energies  are  negative, 
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TAB'LE  6  (continued) 


X^    e-^^^  2,93392  1,91880  1.41779  1.10141  0,64363 

7oAE(i)  02.25  00.13  00.21  00.00  19.55 

fc*-^^  2.94321  1.96909  1.42581  1,11437  0,73923 

7oAE(i)  25,52  60.94  12,98  24,15  77.67 


c  1  \.-     (E.)     2.97295    2.00140    1,48046    1.15508    0,77597 
Solution    i 
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orthogonal  complement  is  insignificant  for  this  state  compared  to  the 
mixing  of  the  ground  state  unperturbed  solution,  i.e.,  the  IIFR  result, 
and  the  perturbation  functions  derived  from  it. 

Since  basis  II  is  very  prejuduced  toward  the  ground  state,  of  the 
four  excited  states,  only  the  first,  the  A  E  state,  behaves  correctly 

at  large  internuclear  separation,  where  it  converges  reasonably  well  to  a 

+  *  1  + 

He  (Is)  +  H(ls)  separated  atom  limit.   The  relationship  between  the  X  E  and 

A  E  curves  is  illustrated  in  Figure  3.   Although  it  is  not  apparent  in 
the  figure,  from  the  full  CI  results  listed  in  Table  7,  the  first  excited 
state  actually  shows  a  slight  minimum  in  the  vicinity  of  R  =  6.5  a.u. 
This  feature  has  been  previously  observed  by  Michels  [66]  who  more  cor- 
rectly positions  the  minimum  at  5.65  a.u. 

In  Figure  4  the  first  and  second  order  reduction  energies  are  com- 
pared with  the  unperturbed  result  and  the  full  CI  solution.   As  before, 

the  extent  of  reduction  given  by  the  low  order  solutions  is  very  satis- 

(2) 

factory.   In  fact,  £    (R)  apparently  converges  to  the  same  separated 

atom  limit  as  the  full  CI,  although  £    (R)  does  not  seem  to  be  able  to 
do  so.   The  eighth  order  result  is  observed  to  be  essentially  identical  to 
the  full  solution. 

As  has  been  mentioned,  the  approximate  quantities  defined  in  (5.32) 
are  obtained  as  an  intermediate  in  the  computational  procedure.   In  order 
for  these  conceptually  appealing  results  to  yield  reasonable  estimates, 
the  interactions  arising  from  the  remaining  p  unperturbed  solutions  need 
to  be  small  as  was  indicated  from  the  ground  state  discussion.   If  so, 
neglecting  the  mixing  should  not  appreciably  affect  the  spectrum  of  the 
results.   That  this  is  actually  the  case  in  this  example  is  illustrated  in 


,.+ 


"A  correlation  diagram  for  the  lower  states  of  Hell  has  been  given 
by  Michels  [66]. 
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Table  8,  where  the  approximate  energies,  {£.   },  and  the  rigorous  selec- 
tively  reduced  energies,  {£.   },  are  compared.   Besides  the  ground  state, 
there  is  no  bounding  property  for  these  approximate  energies,  but  as  is 
observed  from  Table  8,  the  removal  of  the  remaining  reference  functions 
from  the  space  seems  to  have  constrained  the  problem  in  such  a  way  that 
an  effect  very  similar  to  that  perpetrated  by  the  variational  principle 
is  present.   Tlie  resultant  eigenvectors  also  support  the  validity  of  tXils 
approximation. 

An  attempt  was  made  to  obtain  more  realistic  excited  states  for 

J. 
HeH   that  tend  to  separate  correctly  at  infinite  R.   At  a  sacrifice  to 

the  ground  state,  a  third  basis  (III)  was  assumed  that  is  primarily  com- 
posed of  "best"  atom  STO's.   The  basis  is  ls„  (2,00),  2s„  (0.75),  2p   ,, 
^  He  He  o  He 

(0.75),  Is^(l.OO),  2Sj^(0.50),  2p^  ^^(0.50),  ls'^^(3.50),  3d^  He^°'^°^' 

3d^  ^(0.333),  2p^^  He(°-^^)'  ^p;^  He^^'^^^'  ^p^l  H^^'^^^'  ^p^  ^(1.00). 
Since  this  set  contains  functions  which  are  able  to  well  approximate  the 
pertinent  excited  state  separated  atom  limits.  He  (Is),  He(ls2s),  H(ls), 
H(2s),  and  H(2p  )  [66],  the  65  configuration  full  CI  results  at  12  a.u. 
listed  in  Table  9  show  adequate  convergence  to  the  experimental  separated 
atom  values,  -2.500,  -2.1460,  -2.1250,  and  -2.1250  a.u.,  respectively. 
These  curves  are  illustrated  in  Figure  5.   The  change  of  scale  makes  the 


m 


inimum  in  the  A  X!  curve  apparent,  but  the  relatively  broad  and  deep 


1  4- 
minimum  in  the  B  E  curve  is  not  readily  observed  due  to  R  being  limited 

to  8.00  a.u.   From  the  values  in  Table  9  for  larger  R,  though,  this  is 
seen  to  be  the  case.   The  other  two  states  are  repulsive. 

By  comparing  these  calculations  vjith  Michels'  values  [66],  the  depths 
of  the  wells  in  the  two  bound  excited  states  are  still  somewhat  insuffi- 
ciently described  with  basis  III,  though  these  results  are  adequate 


TABLE  8 

Comparison  of  Rigorous  and  Approximate  Multiple  State  Reduction 
Energies  for  ]\e\l^   as  a  Function  of  Internuclear  Separation 

(Ten  function  partitioning;  HE  and  single  excitations;  Basis  II) 
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10.00 
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^(1) 
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2.88522 

2.95820 

2.93533 

2.88739 

2. 8? 568 

0        ^(8) 
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2.88529 
2.89813 
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2.97278 

2.93545 
2.95006 

2.88750 
2.89820 
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2.89733 
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2.89741 
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r(i) 


S' 

1,26363 

1,44918 

1.46924 

1.56815 

1.50428 

^(1) 

2   ^2 

1,26356 

1.44909 

1.46916 

1.56812 

1 , 50420 

^f) 

1,28573 

1.47967 

1.50847 

1.60296 

1.53536 

ef) 

1.28582 

1.48031 

1,51022 

1.60302 

1.53532 

7(1) 

^3 

0.74810 

1.12218 

1,25367 

1.38074 

1,46250 

3   4^' 

0.74800 

1.12206 

1.25356 

1.38062 

1,46244 

ef) 

0.77500 

1.15608 

1.29318 

1.41107 

1,48318 

ef 

0.77485 

1.15461 

1.28917 

1,41099 

1,48317 

TABLE  8  (continued) 
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estimates.   These  excited  states  are  very   atomic-like  and  show  little 

correlation.   This  feature  is  indicated  by  the  CI  eigenvectors  which  show 

appreciable  mixing  for  only  the  first  10  configurations,  all  of  which  are 

obtained  from  the  Is   (2.00)  and  the  nine  possible  sigma  orbitals  which 

He 

may  be  substituted.   Consequently,  to  improve  the  results,  more  extensive 
experimentation  with  the  basis  set  is  warranted  to  give  CI  expansions 
where  the  remaining  configurations  are  of  more  significance. 

A  partitioning  into  the  first  five  configurations  obtained  from  the 
Schmidt  orthogonalized  spin  orbitals  that  arise  from  basis  III  was  made 
to  provide  an  unperturbed  reference  space  for  the  reduction  process.   The 
first  five  elements  of  basis  III  are  the  spatial  orbitals  involved  in 
this  partitioning,  so  the  unperturbed  solution  is  close  to  a  best  atom 
approximation. 

The  reduced  energies  are  not  as  dramatic  in  first  order  as  previously, 

which  probably  signifies  the  lack  of  basis  adaptability.   The  eighth  order 

1  + 
result  is  still  rather  good,  although  it  is  superior  for  the  lowest  A  Z  and 

1  +  1  + 

B  E  states.   The  first,  second,  and  sixth  order  solutions  for  the  B  E 

state  are  illustrated  in  Figure  6.   This  figure  shows  better  agreement  at 

smaller  R,  while  seeming  to  approach  a  somewhat  higher  separated  atom  limit 

than  does  the  full  CI  result. 

This  study  of  the  reduced  partitioning  procedure  has  in  some  ways  been 

preliminary,  as  only  two-electron  molecules  have  been  considered.   Unlike 

many  perturbation  approaches,  though,  the  development  of  tliis  theory  always 

had  as  its  foremost  objective  improved,  or  superior,  treatments  for  larger 

molecules.   As  such,  the  formulation  is  by  no  means  limited  to  very  small 

systems,  but  is  computationally  feasible  for  more  general  cases.   With  the 

success  of  this  investigation,  an  extension  to  more  chemically  interesting 

species  is  justified. 
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